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Abstract 

We consider a communication method, where the sender encodes n classical bits into 1 qubit 
and sends it to the receiver who performs a certain measurement depending on which of the initial 
bits must be recovered. This procedure is called n i-^- 1 quantum random access code (QRAC) 
where p > 1/2 is its success probability. It is known that 2 3^ 1 and 3 1 QRACs (with no 
classical counterparts) exist and that 4 i-^ 1 QRAC with p> 1/2 is not possible. 

We extend this model with shared randomness (SR) that is accessible to both parties. Then 
n 1 QRAC with SR and p > 1/2 exists for any n > 1. We give an upper bound on its success 
probability (the known 2 i^-^ 1 and 3 1 QRACs match this upper bound). We discuss some 
particular constructions for several small values of n. 

We also study the classical counterpart of this model where n bits are encoded into 1 bit instead 
of 1 qubit and SR is used. We give an optimal construction for such codes and find their success 
probability exactly — it is less than in the quantum case. 
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1 Introduction 



1.1 Random access codes 

In general random access code (or simply RAC) stands for "encoding a long message into fewer bits 
with the ability to recover (decode) any one of the initial bits (with some probability of success)" . A 
random access code can be characterized by the symbol "n m" meaning that n bits are encoded into 
m and any one of the initial bits can be recovered with probability at least p. We require that p > 1/2 
since p — 1/2 can be achieved by guessing. In this paper we consider only the case when m = 1. So we 
have the following problem: 

Problem (Classical). There are two parties — Alice and Bob. Alice is asked to encode some classical 
n-bit string into 1 bit and send this bit to Bob. We want Bob to be able to recover any one of the n 
initial bits with high success probability. 

Note that Alice does not know in advance which bit Bob will need to recover, so she cannot send 
only that bit. If they share a quantum channel then we have the quantum version of the previous 
problem: 

Problem (Quantum). Alice must encode her classical n-bit message into 1 qubit (quantum bit) and 
send it to Bob. He performs some measurement on the received qubit to extract the required bit (the 
measurement that is used depends on which bit is needed). 

Both problems look similar, however the quantum version has an important feature. In the classical 
case the fact that Bob can recover any one of the initial bits implies that he can actually recover all of 
them — each with high probability of success. Surprisingly in the quantum case this is not true, because 
after the first measurement the state of the qubit will be disturbed and further attempts to extract 
more information can fail. 

1.2 History and applications 

As noted in [6, 8], the idea behind quantum random access codes or QRACs is very old (relative to 
quantum information standards). It first appeared in a paper by Stephen Wiesner [1] published in 1983 
and was called conjugate coding. Later these codes were re-discovered by Ambainis et al. in [2, 3]. They 

show that there exists 2 1 QRAC and mention its immediate generalization to 3 1 QRAC due 
to Chuang (see also [5] and [8] for more details). However, Hayashi et al. [5] show that it is impossible 
to construct a 4 A 1 QRAC with p > 1/2. We will discuss these results more in Sect. 3.3. 

There has also been work on n m codes with m > 1, see [2, 3, 4]. Ambainis et al. [2] show that if 
a n i-^- TO QRAC with p > 1/2 exists, then to — f2(n/logn), which was later improved by Nayak [4, 3] 
to TO > (1 — H(jp))n, where H(p) = —plogp — (1 — p) log(l — p) is the binary entropy function. Other 
generalizations include: considering c?-valued bits instead of qubits [6, 8] and recovering several rather 
than a single bit [16]. 

Originally quantum random access codes were studied in the context of quantum finite automata 
[2, 3, 4]. However, they also have applications in quantum communication complexity [6, 9, 10, 11], in 
particular for network coding [5, 12] and locally decodable codes [13, 14, 15, 16]. Recently results on 
quantum random access codes have been applied for quantum state learning [17]. 

Experimental feasibility of QRACs and their relation to contextuality and non-locality has been 
discussed in [6, Chapter 7]. Recently a similar protocol called parity- oblivious multiplexing has been 
considered in [7]. It has an additional cryptographic constraint that Alice is not allowed to transmit 
any information about the parity of the input string. In addition [7] also discuss the first experimental 
demonstration of 2 1 and 3 i— > 1 QRACs. 

We want to emphasize the setting in which the impossibility of 4 A- 1 QRAC with p > 1/2 was proved 
in [5]: Alice is allowed to perform a locally randomized encoding of the given string into a one-qubit 
state and Bob is allowed to perform different positive operator-valued measure (POVM) measurements 
to recover different bits. This is the most general setting when information is encoded into a one-qubit 
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state and both parties are allowed to use randomized strategies, but only have access to local coins. 
However, we can consider an even more general setting — when both parties share a common coin. This 
means that Alice and Bob are allowed to cooperate by using some shared source of randomness to agree 
on which strategy to use. We will refer to this source as a shared random string or shared randomness 
(SR). Note that shared randomness is a more powerful resource than local randomness, since parts of 
the shared random string can be exclusively used only by Alice or Bob to simulate local coins. It turns 
out that in this new setting 4 1 QRAC is possible with p > 1/2. In fact, n 1 QRACs with p > 1/2 
can be constructed for all n > 1 (see Sect. 3.8). 

1.3 Outline of results 

In Sect. 2 we study classical n 1 random access codes with shared randomness. In Sect. 2.2 we 
introduce Yao's principle that is useful for understanding both classical and quantum codes. A classical 
code that is optimal for all n is presented in Sect. 2.3.1 and the asymptotic behavior of its success 
probability is considered in Sect. 2.3.2. 

In Sect. 3 we study quantum random access codes with shared randomness. In Sect. 3.3 we discuss 
what is known in the case when shared randomness is not allowed, i.e., 2 i-^ 1 and 3 i-^ 1 QRACs and 
the impossibility of 4 i-^ 1 QRAC. In Sect. 3.6 we give an upper bound of success probabihty of QRACs 
with SR and generalize it in Sect. 3.7 for POVM measurements. In Sect. 3.8 we give two constructions 
of n 1^ 1 QRAC with SR and p > 1/2 for all n > 2 that provide a lower bound for success probability. 

In Sect. 4 we try to find optimal QRACs with SR for several small values of n. In particular, in 
Sect. 4.1 we discuss QRACs obtained by numerical optimization, and in Sect. 4.2 we consider symmetric 
constructions. 

Finally, we conclude in Sect. 5 with a summary of the obtained results (Sect. 5.1), a list of open 
problems (Sect. 5.2) and possible generalizations (Sect. 5.3). 

2 Classical random access codes 

2.1 Types of classical encoding-decoding strategies 

As a synonym for random access code we will use the term strategy to refer to the joint encoding- 
decoding scheme used by Alice and Bob. Two measures of how good the strategy is will be used: the 
worst case success probability and the average success probability. Both probabilities must be calculated 
over all possible pairs (cc, i) where x g {0, 1}" is the input and i G {1, . . . , n} indicates which bit must 
be recovered. We are interested in the worst case success probability, but in our case according to Yao's 
principle (introduced in Sect. 2.2) the average success probability can be used to estimate it. 

Depending on the computational model considered, dilTerent types of strategies are allowed. The 
simplest type corresponds to Alice and Bob acting deterministically and independently. 

Definition. A pure classical n i— > 1 encoding-decoding strategy is an ordered tuple {E, Di, . . . , _D„) that 
consists of an encoding function E : {0, 1}" i— > {0, 1} and n decoding functions Di : {0, 1} i— > {0, 1}. 

These limited strategies yield RACs with poor performance. This is because Bob can recover all bits 
correctly for no more than two input strings, since he receives either or 1 and acts deterministically 
in each case. For all other strings at least one bit will definitely be recovered incorrectly, therefore 
the worst case success probability is 0. If we allow Alice and Bob to act probabilistically but without 
cooperation, then we get mixed strategies. 

Definition. A mixed classical n i—f I encoding-decoding strategy is an ordered tuple {Pe, Pon ■ ■ ■ , Pd„) 
of probability distributions. Pe is a distribution over encoding functions and P/j. over decoding func- 
tions. 

It is obvious that in this setting the worst case probability is at least 1/2. This is obtained by 
guessing — we output either or 1 with probability 1/2 regardless of the input. Formally this means 
that for each i, Fq. is a uniform distribution over two constant decoding functions and 1. It has been 
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shown that in this setting for 2 i— s- 1 case one cannot do better than guessing, i.e., there is no 2 i— > 1 
RAC with worst case success probabihty p > 1/2 [2, 3]. 

However, we can allow cooperation between Alice and Bob — they can use a shared random string 
to agree on some joint strategy. 

Definition. A classical n i—s- 1 encoding-decoding strategy with shared randomness is a probability 
distribution over pure classical strategies. 

Note that this is the most general randomized setting, since both randomized cooperation and local 
randomization are possible. This is demonstrated in the following example. 

Example. Consider the following strategy: randomly agree on i e {1, . . . ,n} and send the ith bit; if 
the ith bit is requested, output the received bit, otherwise guess. This strategy can formally be specified 
as follows: uniformly choose a pure strategy from the set 

IJ {(e^jCi, . . .,Ci-i,d,gi, . . ..gn-i) \ c G {do,di}'"\g G {do, (1) 

iG{l,...,n} 

where the encoding function is given by ei{x) = Xi and decoding functions do, ^i, and d are given by 
doih) = 0, di(h) = 1, and d(h) = 6, where h is the received bit. The total amount of required randomness 
is n — 1 + log n bits, because one out of n ■ 2^~^ pure strategies must be selected. Note that only log n of 
these bits must be shared among Alice and Bob, so that they can agree on the value of i. The remaining 
n — 1 random bits are needed only by Bob for choosing random decoding functions c e {do, di}' ^ and 
g e {do,di}""*. 

Note that the amount of randomness used in the above example can be reduced. Since only one bit 
must be recovered, there is no need to choose each of the decoding functions independently. Thus Bob 
needs only one random bit that he will output whenever some bit other than the ith bit is requested. 
This is illustrated in the next example. 

Example. Alice and Bob uniformly sample a pure strategy from the following set: 

{(ei,c^^^^,d,c, . . .^) I 1 < i < 7i,c e {do,di}}. (2) 

i—l n—i 

This requires log n random bits to be shared among Alice and Bob and 1 private random bit for Bob, 
i.e., 1 + logn random bits in total. 

We are interested in classical strategies with SR, because they provide a classical analogue of QRACs 
with SR. However, in this setting finding the optimal strategy seems to be hard, therefore we will turn 
to Yao's principle for help. 

2.2 Yao's principle 

When dealing with randomized algorithms, it is hard to draw general conclusions (like proving optimality 
of a certain randomized algorithm) because the possible algorithms may form a continuum. In such 
situations it is very helpful to apply Yao's principle [18]. This allows us to shift the randomness in the 
algorithm to the input and consider only deterministic algorithms. 

Let iS be a classical strategy with SR. One can think of it as a stochastic process consisting of applying 
the encoding map E to the input x, followed by applying the decoding map Di to the ith bit. Both of 
these maps depend on the value of the shared random string. The result of S is S{x,i) = Di{E{x)), 
which is a stochastic variable over the set {0, 1}. Let PT[S{x,i) — Xi] denote the probability that the 
stochastic variable S{x, i) takes value Xi. Then the worst case success probability of the optimal classical 
strategy with SR is given by 

maxminPr[5(a;, i) = a;i]. (3) 
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Let fi be some distribution over the input set {0, 1}" x {1, . . . , n} and let Pr^[P(a;, i) — Xi] denote the 
expected success probabihty of a pure (deterministic) strategy V. If the "hardest" input distribution is 
chosen as /x, then the expected success probability of the best pure strategy for this distribution is 

minmaxPr^[7'(a;, i) = Xi]. (4) 

Yao's principle states that the quantities given in (3) and (4) are equal [18]: 

maxminPr[5(x, i) = Xi] = minmaxPru[7'(a:, «) = Xi\. (5) 

Thus Yao's principle provides us with an upper bound for the worst case probability (3). All we have 
to do is to choose an arbitrary input distribution /ip and find the best pure strategy Vq for it. Then 
according to Yao's principle we have 

PVo[''^o(a;, *) = Xi] > maxminPr[5(a;,i) = x.^], (6) 

with equality if and only if fj,Q is the "hardest" distribution. It turns out that for random access codes 
the uniform distribution rj is the "hardest". To prove it, we must first consider the randomization 
lemma. 

Lemma 1. yV3S : min^^^i Pr[5(a;, i) = Xi] = Pr^[7'(x,i) = Xi], where rj is the uniform distribution. In 
other words: the worst case success probability of S is the same as the average case success probability 
of V with uniformly distributed input. 

Proof. This can be achieved by randomizing the input with the help of the shared random string. 
Alice's input can be randomized by XOR-ing it with an n-bit random string r. But Bob's input can 
be randomized by adding (modulo n) a random number d € {0, . . . , n — 1} to it (assume for now that 
bits are numbered from to n — 1). To obtain a consistent strategy, these actions must be identically 
performed on both sides, thus a shared random string of n + logn bits^ is required. Assume that E 
and Di are the encoding and decoding functions of the pure strategy V] then the new strategy S is 

E'{x)=E{S\iiiU{x®r)), (7) 

mod 71 

{b)®r„ (8) 

where Shift(i(s) substitutes Si+rfmod n by Si in string s. Due to input randomization, this strategy has 
the same success probability for all inputs {x,i), namely 

n — 1 ^ 

Pr[Six,i)=x,] = E^;rT7^P^-[^(y'J') = %]=P^-''[^(2^'^')=?/j]' (9) 

yG{0,l}" j=0 

coinciding with the average success probability of the pure strategy V. □ 

Now we will show that inequality (6) becomes an equality when /iq = Vi meaning that the uniform 
distribution rj is the "hardest" . 

Lemma 2. The minimum of (4) is reached at the uniform distribution rj, i.e., 

min max Pr^ [7^(2;, i) — Xi] — maxPr^[7'(a;, i) — Xi\. (10) 

Proof. From the previous Lemma we know that there exists a strategy with SR Sq such that 

minPr[5o(a;, i) = Xi] = max Pr„ [7^(2;, i) ~ Xi] (11) 



^We will not worry about how Bob obtains a uniformly distributed d from a string of random bits when n ^ 2^ 
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(5o is obtained from the best pure strategy by prepending it with input randomization). However, 
among all strategies with SR there might be one that is better than Sq, thus 

mjixminPr[5(x, i) = Xi] > maxPr.,j7'(a;, i) = Xi]. (12) 

But if we put — 7] into inequality (6), we obtain 

max Pr„ [7^(2;, i) = Xi] > maxminPr[iS(a;, i) = Xi], (13) 

V S x.i 

which is the same as (12), but with reversed sign. This means that both sides are actually equal: 

maxPr„[7^(a;, i) = Xj] — maxminPr[iS(a;, i) — Xi\. (14) 

V S x,i 

Applying Yao's principle to the right hand side of (14) we obtain the desired equation (10). □ 
Theorem 1. For any pure strategy V 

Frjj[V{x,i) — Xi] < maxminPr[5(a;, i) — Xi], (15) 

5 x,i 

with equality if and only if V is optimal for the uniform distribution rj. 

Proof. To obtain the required inequality, do not maximize the left hand side of equation (14), but put 
an arbitrary V. It is obvious that we will obtain equality if and only if V is optimal. □ 

This theorem has important consequences — it allows us to consider pure strategies with uniformly 
distributed input rather than strategies with SR. If we manage to find the optimal pure strategy, then 
we can also construct an optimal strategy with SR using input randomization^. If the pure strategy is 
not optimal, then we get a lower bound for the strategy with SR. 



2.3 Classical n 1 RAC 

Before considering n 1-^ 1 QRACs with shared randomness, we will find an optimal classical n i— > 1 
RAC with shared randomness and derive bounds for it. 



2.3.1 Optimal strategy 

According to Theorem 1 we can consider only pure strategies. As a pure strategy is deterministic, for 
each input it gives either a correct or a wrong answer. To maximize the average success probability we 
must find a pure strategy that gives the correct answer for as many of the n ■ 2" inputs as possible — such 
a strategy we will call an optim,al pure strategy. 

Let us first consider the problem of finding an optimal decoding strategy, when the encoding strategy 
is fixed. An encoding function E : {0, 1}" 1^ {0, 1} divides the set of all strings into two parts: 

Xo^{xe{0,ir\E{x) = 0}, 

X, = {xe{0,ir\E{x)^l}. ^ ' 

If Bob receives bit b, he knows that the initial string was definitely from the set Xi,, but there is no way 
for him to tell exactly which string it was. However, if he must recover only the ith bit, he can check 
whether there are more zeros or ones among the ith bits of strings from set Xh. More formally, we can 
introduce the symbol N^{k) that denotes the number of strings from set Xi, that have the bit k in ith 
position: 

N^{k)^\{xeXb\x,^k}\, (17) 

^If the encoding function depends only on the Hamming weight of the input string x (e.g., majority function) and the 
decoding function does not depend on i, there is no need to randomize over j, so n instead of n + logn shared random 
bits are enough. 
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Therefore the optimal decoding strategy Di : {0,1} {0,1} for the ith bit is 



if N^(0) > N''(l), 
D,{b) ^{ ' - (18) 

I 1 otherwise. 

Of course, if N^{0) — N^{1), Bob can output 1 as well. For pure strategies there are only 4 possible 
decoding functions for each bit: 0, 1, b, or NOT 6. But this is still quite a lot so we will consider the 
two following lemmas. The first lemma will rule out the constant decoding functions and 1. 

Lemma 3. For any n there exists an optimal pure classical n i— > 1 RAC that does not use constant 
decoding functions and 1 for any bits. 

Proof. We will show that if there exists an optimal strategy that contains constant decoding functions 
for some bits, then there also exists an optimal strategy that does not. Let us assume that there is an 
optimal strategy with constant decoding function for the ith bit (the same argument goes through 
for 1 as well). Then according to equation (18) we have 7V°(0) > Nf{l) and 7V/(0) > Nl{l). Note that 
A^°(0) + A^/(0) = 7V°(1) + 7V/(1) = 2"-\ because Xi=0\n exactly half of all 2" strings. This means 
that actually iV°(0) = iV°(l) and iV/(0) = iV/(l). If we take a look at (18) again, we see that in such 
situation any decoding strategy is optimal and we can use any non-constant strategy instead. □ 

Lemma 4. For any n there exists an optimal pure classical n t-^ 1 RAC that does not use decoding 
function NOT b for any bits. 

Proof. We will show that for each pure strategy V that uses negation as the decoding function for the 
ith bit, there exists a pure strategy V' with the same average case success probability that does not. 
If V consists of encoding function E and decoding functions Dj, then V' can be obtained from V by 
inverting the zth bit before encoding and after decoding: 

E'{x) = E(NOT^x), (19) 

m - i''''^.'''^'^ (20) 

\Dj[b) otherwise, 

where NOT^ inverts the ith bit of string. It is obvious that V and V' have the same average success 
probabilities, because if V gives the correct answer for input (x, i) then V' gives the correct answer for 
input (NOTi x,i). The same holds for wrong answers. □ 

Theorem 2. The pure classical n t—i- 1 RAC with identity decoding functions and majority encoding 
function is optimal. 

Proof. According to Lemma 3 and Lemma 4, there exists an optimal pure classical n i— > 1 RAC with 
identity decoding function for all bits. Now we must consider the other part — finding an optimal 
encoding given a particular (identity) decoding function. It is obvious that in our case optimal encoding 
must return the majority of bits: 

E'ix)=l' ^^l^l<"/2' (21) 
I 1 otherwise, 

where is the Hamming weight of string x (the number of ones in it). □ 
2.3.2 Asymptotic bounds 

Let us find the exact value of the average success probability for the optimal pure RAC suggested in 
Theorem 2. We will separately consider the even and odd cases. 

In the odd case (n = 2m + 1) the average success probability is given by 




^(^-+^)= (2m+l).2^»-+i '^ >^ ^ 11' (22) 
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p(n) 



10 11 12 13 14 15 




Figure 1: Exact probability of success p{n) for optimal pure classical n i— > 1 RAC (black dots) according 
to (24) and its approximate value (dashed line) according to (27). Dotted lines show upper and lower 
bounds of p{n) for odd and even n according to inequalities (29) and (30). 



where the factor 2 stands for either zeros or ones being the majority, and ('^'";^^) stands for the number 
of strings where the given symbol dominates and appears exactly i times. 

In the even case (n = 2m) there are a lot of strings with the same number of zeros and ones. These 
strings are bad, because with majority encoding and identity decoding it is not possible to give the 
correct answer for more than half of all bits. The corresponding average success probability is given by 



Pi^rn)^—— 2 V ^ ^ +m[ , (23) 




where the last term stands for the bad strings. 

In Appendix A we give a combinatorial interpretation of the sums in (22) and (23). Equations (143) 
and (144) derived in Appendix A can be used to simplify p(2m + 1) and p(2m), respectively. It turns 
out that both probabilities are equal: 

1 1 ^2m\ 

p{2m) = p{2m +1) = -^ + ^ J • (24) 

These two expressions can be combined as follows: 

We can apply Stirling's approximation [20] to! ~ iJ^Y^ \J2-Km to (24) and obtain 

p(2m) = p{2ra + 1) « ^ + (26) 

2 2y TTTO 

If we put TO w |, then (26) turns to 

p(n)«J + -i=. (27) 
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We see that the value of (27) approaches 1/2 as n increases. Thus the obtained codes are not very 
good for large n, since p = 1/2 can be obtained by guessing. We will observe a similar (but slightly 
better) behavior also in the quantum case. The exact probability (24) and its approximation (27) are 
shown in Fig. 1. 

For odd and even cases asymptotic upper and lower bounds on p{n) can be obtained using the 
following inequaUty [20]: 



/ / Tl\^ 1 ; / Tl\^^ 1 

V2Tm\^-j eT2M^ < n! < V27rn (^-j . (28) 
For the odd case we have 



y 12ri-ll 



exp I T2;r^ - 6;r36 _ i exp 



, <p(n)--< ^ , ^, (29) 

v/2^(n - 1) 2 ^27r(n-l) 



but for the even case 



exp 



(l2n 6n+l) 1 ^'''■P ( 12n+l 



\/2TTn 2 v27rn 

All four bounds are shown in Fig. 1. 



<pH-o< -■ (30) 



3 Quantum random access codes 
3.1 Visualizing a qubit 

When dealing with quantum random access codes (at least in the qubit case), it is a good idea to try 
to visualize them. We provide two ways. 

3.1.1 Bloch sphere representation 

A pure qubit state is a column vector £ C^. It can be expressed as a linear combination lip) — 
a\0) +/3|1), where |0) = (J) and |1) = (o). The coefficients a, (3 S C must satisfy \a\^ + \(3\^ = 1. 
Since the physical state is not affected by the phase factor (i.e., {tp) and e*"^ are the same states for 
any g M) , without the loss of generality one can write 

where < 6 < tt and < ip < 27t (the factor 1/2 for 9 in (31) is chosen so that these ranges resemble 
the ones for spherical coordinates in M.^). 

For almost all states jV') there is a unique way to assign the parameters and (p. The only exceptions 
are states |0) and |1), that correspond to 9 = and 9 = n, respectively. In both cases (p does not affect 
the physical state. Note that the spherical coordinates with latitude 9 and longitude ip have the same 
property, namely — the longitude is not defined at poles. This suggests that the state space of a single 
qubit is topologically a sphere. 

Indeed, there is a one-to-one correspondence between pure qubit states and the points on a unit 
sphere in M'^. This is called the Bloch sphere representation of a qubit state. The Bloch vector for state 
(31) is r = {x, y, z), where the coordinates (see Fig. 2) are given by 

X = sin cos (p, 

y = sin 6* sin (32) 
z — cos 9. 

Given the Bloch vector r — {x, y, z), the coefficients of the corresponding state — a \ 0) + (3 \1) can 
be found as follows [19, pp. 102]: 
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Figure 2: Angles 9 and (p of the Bloch Figure 3: Geometric interpretation of or- 

vector corresponding to state 1-0). thogonal measurement. 



with the convention that (0,0, —1) corresponds to a = and (3=1. 

The density matrix of a pure state {ip) is defined as p = For the state \ip) in (31) we have 

1/1+ cos 61 e-''^sme\ 1,, 
P=2[e^^sine l-cos0) = 2^^ + ''^^^y''^ + '^'^' ^^^^ 

where (a;, y, z) are the coordinates of the Bloch vector r given in (32) and 

^=(J J), '^x=(J J), ^.= (- V)' ^^=(o -i) 

are called Pauli matrices. We can write (34) more concisely as 

P=\{I + r-a) (36) 

where r = {x,y,z) and cr = (cr^r, cr^, ct^). 

If ri and r2 are the Bloch vectors of two pure states \ipi) and IV'2), then 

|(Vi|^2)|' - Tr(pip2) = ^(1 + ri • ra). (37) 

This relates the inner product in C'^ to the one in M.^. Since ri and r2 are unit vectors, Ti • r2 = cos a, 
where a is the angle between ri and r2. 

An orthogonal measurement M on a qubit can be specified by a set of two orthonormal states: 
M = {IV'o) I l^i)}- Orthonormality means that {^pilijjj) = 6ij. If we measure a qubit that is in state {ip) 
with measurement M then the outcome will be either or 1 and the state will "collapse" to jV'o) or 1-01) 
with probabihties KV'olV') I and |(V'i|'0)l i respectively. Observe that for orthogonal states equation (37) 
implies ri ■ r2 = —1, therefore they correspond to antipodal points on the Bloch sphere. If we denote 
the angle between the Bloch vectors of \ip) and \ipo) by a, then according to (37) the probabilities of 
the outcomes are ^ 

Po = ;^(1 + cosa), 

I (38) 
Pi = - cosa). 

There is a nice geometrical interpretation of these probabilities. If we project the Bloch vector corre- 
sponding to IV') on the axes spanned by the Bloch vectors of \tpo) and l^i) (see Fig. 3), then po = di/2 
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Figure 4: Curves of constant 9 and ip before (on the left) and after the Hadamard transformation (on 
the right). Initially the curves of constant 9 are concentric circles, but after the transformation they 
appear as deformed circles around both poles. The curves of constant ip transform form radial rays to 
"field lines" connecting both poles. The image on the left appears to have only the North pole |0), since 
the Bloch sphere is punctured at the South pole |1) which must be identified with the boundary of the 
unit disk. The "left pole" and "right pole" in the image on the right correspond to the states |1) and 
|0), respectively. 



and pi = do/2 (note the different indices), where do is the distance between the projection and |V'o)j 
but di is the distance between the projection and IV'i)- Observe that vectors on the upper hemisphere 
have greater probability to collapse to \tpo), but on lower hemisphere, to IV'i)- On the equator both 
probabilities are equal to ^. 

3.1.2 Unit disk representation 

There is another way of visualizing a qubit. Unlike the Bloch sphere representation, this way of repre- 
senting a qubit is not known to have appeared elsewhere. The idea is to use only one complex number 
to specify a pure qubit state |-!/)) = ( ^ ) S C^. This is possible since can be written in the form (31), 
which is completely determined by its second component 



identified with a unit disk in the complex plane (the polar coordinates assigned to axe (r, ip), where 
r — sin |). The origin (3 = corresponds to lip) = |0), and all points on the unit circle \(3\ = 1 are 
identified with \tp) = |1), since e**^ |1) corresponds to the same quantum state for all G K. 

The relation between the unit disk representation and the Bloch sphere representation can be visu- 
alized as follows: 

• the unit disk is obtained by puncturing the Bloch sphere at its South pole and flattening it, 

• the Bloch sphere is obtained by gluing together the boundary of the unit disk. 

It is much harder to visualize how a unitary transformation acts in the unit disk representation. Let 
us consider a simple example. 




(39) 



The first component 




As 1/3 1 < 1, the set of all possible qubit states can be 
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Example. Let us consider the action of the Hadamard gate H = ( i -i ) ™ the unit disk represen- 
tation. Note that = I thus H is an involution (self-inverse). It acts on the standard basis states as 
follows: 

^|0) = 7llO) + 7l|l) = l+)' (40) 
^|l> = ;)5|0)-;i3|l) = |-). (41) 

The way H transforms the curves of constant 9 and (p is shown in Fig. 4. From equation (40) we see 
that the origin /3 = corresponding to |0) is mapped to the "right pole" /3 = ^ corresponding to |-|-) 
(and vice versa). Recall that all points on the boundary of the unit disk in Fig. 4 (on the left) are 
identified with |1). Thus equation (41) tells us that the unit circle = 1 is mapped to the "left pole" 
(3 = Fig. 4 (on the right) corresponding to |— ) (and vice versa). This means that |— ) is mapped 

to the boundary of the unit disk in Fig. 4 (on the right). 

Since we use only one complex number /3 to represent a quantum state, a finite set of quantum 
states {/3i, /32, • ■ • , Pn} can be represented by a polynomial 

c(/3-/3i)(/3-/32)---(/3-/3„) (42) 

whose roots are (3i (here c ^ is arbitrary). We will use this representation in Sects. 3.3 and 4.1 to 
describe the qubit states whose Bloch vectors are the vertices of certain polyhedra. It is surprising that 
for those states the values of c can be chosen so that the resulting polynomials have integer coefficients. 



3.2 Types of quantum encoding-decoding strategies 

Let us now consider the quantum analogue of a pure strategy. 

Definition. A pure quantum n i— > 1 encoding-decoding strategy is an ordered tuple (E, Mi, . . . , M„) that 
consists of encoding function i? : {0, 1}" i— > and n orthogonal measurements: Mi = {jV'o) ' 

If Alice encodes the string x with function E, she obtains a pure qubit state 1-0) — E{x). When Bob 
receives jV') and is asked to recover the ith bit of x, he performs the measurement Mi. The probability 
that Bob recovers Xi correctly is equal to 

p(x,i)= KV'^JV)!'- (43) 

As in the classical setting, we can allow Alice and Bob to have probabilistic quantum strategies 
without cooperation. Though we will not need it, mixed quantum strategies can be defined in complete 
analogy with mixed classical strategies. 

Definition. A mixed quantum 7i i— > 1 encoding- decoding strategy is an ordered tuple (-Pb, Pmi t ■ ■ i Pm„) 
of probability distributions. Pe is a distribution over encoding functions E and are probability 
distributions over orthogonal measurements of qubit. 

The main objects of our research are quantum strategies with cooperation, i.e., with shared ran- 
domness. They are defined in complete analogy with the classical ones. 

Definition. A quantum n I encoding-decoding strategy with shared randomness is a probability 
distribution over pure quantum strategies. 

We would like to point out two very important things about quantum strategies with shared ran- 
domness. The first thing is that all statements about classical strategies with SR in Sect. 2.2 are valid 
for quantum strategies as well (the only difference is that "pure strategy" now means "pure quantum 
strategy" instead of "pure classical strategy" and "strategy with SR" means "quantum strategy with 
SR" instead of "classical strategy with SR"). The most important consequence of this observation is 
that Theorem 1 is valid also for quantum strategies with SR. This means that the same technique of 
obtaining the upper bound can be used in the quantum case, i.e., we can consider the average success 
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Figure 5: Bloch sphere representation of Figure 6: Bloch sphere representation of 

encoding for 2 i-^ 1 quantum random ac- encoding for 3 i-^ 1 quantum random ac- 

cess code. cess code. 



probabihty of a pure quantum strategy instead of the worst case success probabihty of the quantum 
strategy with SR. 

The second thing is that the quantum strategy with SR is the most powerful quantum encoding- 
decoding strategy when both kinds of classical randomness (local and shared) is allowed. However, it is 
not the most general strategy, since it cannot be used to simulate certain classical strategies, e.g., the 
ones with fixed output. However, it turns out that the ability to simulate such strategies does not give 
any advantage (see Sect. 3.7 and Appendix B). 



3.3 Known quantum RACs 

In [2, 3] it has been shown that for 2 i-^ 1 classical RACs in the mixed setting the decoding party 
cannot do better than guessing, i.e., the worst case success probability cannot exceed 1/2. However, 
if quantum states can be transmitted, there are pure quantum 2 i— > 1 and 3 i-^- 1 schemes [2, 3]. This 
clearly indicates the advantages of quantum RACs. On the other hand, a quantum 4 i— s- 1 scheme 
cannot exist [5]. We will review these results in the next three sections. 



3.3.1 The 2 1 QRAC 

The 2 1-^- 1 QRAC is described in [2, 3, 5]. The main idea is to use two mutually orthogonal pairs of 
antipodal Bloch vectors for measurement bases. For example, let Mi and M2 be the measurements along 
the X and y axes, respectively. The corresponding Bloch vectors are Vi = (±1, 0, 0) and V2 = (0, ±1,0). 
The measurement bases are 

"-{;^C)';^(-.)}- 

The planes orthogonal to the x and y axes cut the Bloch sphere into four parts. Note that in each part 
only one definite string can be encoded (otherwise the worst case success probability will be less than ^ ) . 
According to (38), all encoding points must be as far from both planes as possible in order to maximize 
the worst case success probability (recall the geometrical interpretation of the measurement shown in 
Fig. 3). In our case the best encoding states are the vertices of a square -^{±1, ±1, 0) inscribed in the 
unit circle on the xy plane (see Fig. 5). Given a string x = X1X2, the Bloch vector of the encoding state 
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can be found as follows: 

1 /(-I)" 

r{x) = — (-1)-^ I . (46) 



V2 

The corresponding encoding function is 







E{x„x,) = ^ |0) + ^ ' ' |1> . (47) 

The success probability is the same for all input strings and all bits to be recovered: 

p= ^ (i + costI = ^ + ^ « 0.853553. (48) 
2 V 4/2 2V2 



3.3.2 The 3 ^ 1 QRAC 



It is not hard to generalize the 2 1 QRAC to a 3 1-^ 1 code — just take three mutually orthogonal pairs 
of antipodal Bloch vectors, i.e., the vertices of an octahedron [5, 8]. The third pair is — (0,0, ±1) 
and the corresponding measurement basis is 



(49) 



In this case we have three orthogonal planes that cut the sphere into eight parts and only one string 
can be encoded into each part. In this case the optimal encoding states correspond to the vertices of a 
cube ^(±1, ±1, ±1) inscribed in the Bloch sphere (see Fig. 6). The Bloch vector of the encoding state 
of string x — X1X2X3 is 



r{x) 




{-ir 



The corresponding encoding function is E{xi,X2, X3) 
given by 



(50) 



' /3 |1) with coefficients a and /3 explicitly 



a 




i-iy 



2^/3 ' 



(51) 



y^6 + 2\/3(-l)^3 

In fact, the coefficients (3 are exactly the eight roots of the polynomial'' 

36/3* + 24/3'* + 1 

This code also has the same success probability in all cases: 

1 1 



P 



2^/3 



0.788675. 



(52) 



(53) 



3.3.3 Impossibility of the 4 1 QRAC 

Hayashi et al. [5] have shown that 2 1^ 1 and 3 1-^ 1 codes discussed above cannot be generalized for 4 
(and hence more) encoded bits. The reason is simple — it is not possible to cut the Bloch sphere into 16 
parts with 4 great circles (see the proof below). Thus the number of strings will exceed the number of 
parts, hence at least two strings must be encoded in the same part. This makes the worst case success 
probability drop below ^. 

^The unit disk representation of a quantum state and the representation of a finite set of quantum states using a 
polynomial was discussed in Sect. 3.1.2. 
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Figure 7: Gnomonic projection transforms great circles to lines and vice versa. 



A 




Figure 8: Cutting the plane with 3 lines 
into 7 parts. 



Figure 9: Cutting the sphere with 4 great 
circles into 14 parts (seven diametrically 
opposite parts are equal). 



Let us consider how many parts can be obtained by cutting a sphere with 4 great circles. Without 
loss of generality we can assume that the first great circle coincides with the equator. We use the 
gnomonic projection (from the center of the sphere) to project the remaining three circles to a plane 
tangent to the South pole. Note that great circles are transformed into lines and vice versa (see Fig. 7), 
thus we will obtain three lines. Also note that each region in the plane corresponds to two (diametrically 
opposite) regions on the sphere. It is simple to verify that three lines cannot cut the plane into more 
than 7 parts (see Fig. 8). Thus the sphere cannot be cut into more than 14 parts with four great 
circles.* An example achieving the upper bound is shown in see Fig. 9 (see also Figs. 25 and 26). Using 
essentially the same argument for generalized Bloch vectors Hayashi et al. [5] show that 2^™ i-^ m 
QRACs with p > 1/2 do not exist for all to > 1. The generalized Bloch vector will be briefly introduced 
in Sect. 5.3. 

3.4 Optimal encoding for given decoding strategy 

We just reviewed the known results on pure n >—>■ I quantum random access codes. From now on we 
will consider QRACs with shared randomness. In this section we will show how to find the optimal 
encoding strategy for a given decoding strategy. More precisely, we will show that the measurement 

''In general, if we have n great circles on the sphere, the maximal number of parts we can obtain is twice what we can 
obtain by cutting the plane with n — 1 lines. If each line we draw intersects all previous lines and no three lines intersect 
at the same point, the sphere is cut into n(n — 1) + 2 parts after the inverse gnomonic projection. 
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directions of a QRAC with SR determine the corresponding optimal encoding states in a simple way. 

An orthogonal measurement for the ith bit is specified by antipodal points on the Bloch sphere: 
Mi — {vi,—Vi}. Let Vx be the Bloch vector that corresponds to the quantum state in which string 
X e {0, 1}" is encoded. According to equations in (38) the success probability for input {x,i) is 

p{x,i) = ^{l + {-ir'v,-rx) (54) 
and the average success probability is given by 



1 " 1 




1 \ (55) 

E r-.-E("ir'-0- 



re{oa}" i=i 



In order to maximize the probability p, we only need to maximize S^^j. in equation (55) over all possible 
measurements Vi and encodings (in total n + 2" unit vectors in M^). We will denote the maximum 
of S^y.r by S(ji): 

n 

S(n) = max Sj, r ^ max } max r^-/ {—lY^Vi. (56) 

{■>'^},{rA ' {■>'^} ^f--;,„ ^ 

a;e{0,l}" i=l 

If we define 

= E(-l)'''^" (57) 

then it is obvious that the scalar product ■ in (56) will be maximized when is chosen along the 
same direction as v^, i.e. Tx = v^j when 7^ 0. In this case we have ■ Vx = \\Vx\\ and 



S{ji) = max 2. 



£ce{o,i}" 



E(-i)''^' 



(58) 



Therefore we only need to maximize over all possible measurements succinctly represented by n unit 
vectors Vi € M"^, because the optimal encoding is already determined by measurements (see Sect. 4.1 
for some numerical results obtained in this way). When the value of S(ri) is found, then according to 
(55) the corresponding probability is 

.(..)^i(i + S). (-) 

We can observe a connection between quantum and classical RACs with SR. Assume that Marge 
and Homer^ have to implement n 1— > 1 QRAC with SR and are deciding what strategies to use — Homer 
is responsible for choosing the measurements, but Marge has to choose how to encode the input string. 
Once they have decided, they have to follow the agreement and cannot cheat. Unfortunately, Homer is 
foolish and he proposes to measure all bits in the same basis. Luckily Marge is clever enough to choose 
the optimal encoding for Homer's measurements. According to the discussion above, she has to use the 
majority encoding function. Thus the obtained QRAC is as good as an optimal classical RAC discussed 
in Sect. 2.3.1, Theorem 2. 

It looks plausible that using the same measurement for all bits is the worst decoding strategy. 
However, we have not proved this, so we leave it as a conjecture: 

Conjecture. For any choice of measurements there is an encoding such that the resulting n i—f 1 
quantum RAC with SR is at least as good as the optimal n 1— > 1 classical one. 



^In this scenario it is more convenient to replace Alice and Bob with Marge and Homer from The Simpsons. 
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3.5 Relation to a random walk in 



QRACs with shared randomness are related to random walks in M^. This relation can be seen by suitably 
interpreting equations (58) and (59). Let us consider an rt i-^ 1 QRAC with SR whose measurement 
directions are given by unit vectors {vi} and let us assume that the corresponding optimal encoding 
for these measurements is used as described in the previous section. Then we can write the success 
probability p{vi, . . . , t)„) of this QRAC in the following suggestive form: 



p{vi 



,■"«) = + ^ d{vi,. . . ,vn)y 



where 



d{vi, . . . , v„) 



E 



1=1 



(60) 



(61) 



is the average distance traveled by a random walk whose ith step is Vi or —Vi, each with probability 
1/2. For example, Vi = V2 = ■ ■ ■ = Vn corresponds to a random walk on a line and d(vi, . . . , «„) is the 
average distance traveled after n steps of this walk. Recall from the previous section that this choice of 
{vi} corresponds to the optimal classical RAC and we conjecture that this is the worst possible choice. 
Similarly, if we choose roughly one third of vectors {vi} along each coordinate axis, we obtain a random 
walk in a cubic lattice and d{vi, . . . , v„) is the average distance traveled when roughly n/3 steps are 
performed along each coordinate axis (see Sect. 3.8.2). 

In Sects. 3.8 we will use this relation between random access codes and random walks to prove a 
lower bound for the success probability of n 1 QRACs with SR. 



3.6 Upper bound 

In this section we will derive an upper bound for S{n). For this purpose we rewrite the equation (58) 
in the following form: 

(62) 



S(n) — max S-^ 



where 



ae{i,-iy- 



i=l 



(for convenience we take the sum over the set {1, —1}" instead of {0, 1}"). 
Lemma 5. For any unit vectors Vi, . . . ,Vn we have 

^ llai^iH 1- a„t;„||^ = n • 2". 

ai,...,a„e{l,-l} 



(63) 



(64) 



l-^ill =2. 



Proof. For n = 1 we have 

E 

ai6{l -1} 

Let us assume that equation (64) holds for n — k. Then for n = fc + 1 we have 



^ llai'UiH 

ai,...,afc,afc + ie{l,-l} 

If we write out the sum over Uk+i explicitly, we obtain 



a-kVk + ak+iVk+i\ 



^ (||ait;i H h flfcf fe + ■Ufc+ill^ + \\aivi H h akVk - Vk+i\ 



(65) 



(66) 



(67) 



ai,...,afce{l,-l} 
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We can use the parallelogram identity 



Ml + M2I 



ll«2| 



which holds for any two vectors Ui and 1x2, to simplify the sum as follows: 



E 2 

ai,...,afc6{l,-l} 



\Vk+l\ 



(68) 



(69) 



We know that v^+i is a unit vector and we have assumed that (64) holds for n = k; therefore (69) 

simplifies to 2 (fc • 2*= + 2*^) = (fc + 1) • 2*^+1. □ 

We will use the previous lemma to obtain an upper bound for defined in (63). According to (62) 
this will give us an upper bound for S{n) as well. 



Lemma 6. For any set of unit vectors {■Ui}"^]^, the inequality S'„ < ^/n ■ 2" holds. 

Proof. We can interpret the first sum in equation (63) as an inner product with (1, , 
the Cauchy-Schwarz inequality x ■ y < \\x\\ \\y\\ says that 



,i)e 



p2" 



\ ae{i-ir 



Oil), 



= \/2" Vn ■ 2" = V^- 2", 



where Lemma 5 was used to obtain the first equality. 

Theorem 3. For any rn^ I QRAC with shared randomness, p < h 



Then 



(70) 



□ 



1 

2yS' 



Proof. From Lemma 6 we have < \/n ■ 2". From equation (62) we see that the same upper bound 
applies to S{n). Putting this into (59) we get 



P< 



□ 



In particular, this means that the known 2 1 and 3 i— > 1 QRACs discussed in Sect. 3.3 cannot be 
improved even if shared randomness is allowed. 

The intuition behind this upper bound is as follows. If instead of M.^ the Bloch vector of a qubit 
state would be in K", we could choose all n measurements to be mutually orthogonal. For example, 
we could choose the vectors forming measurement bases to be the vertices of the cross polytope, i.e., all 
permutations of (±1, 0, . . . , 0). The optimal encoding corresponding to this choice are the vertices of 
the hypercube, i.e., points (±1,±1, . . . ,±1), thus all terms in equation (63) are equal to y/n and sum 
to 2'^y/n, so the probability (59) is |(1 + %^) = | + 2^' have only three dimensions, the 

actual probability should not be larger. 



3.7 General upper bound 

Let us prove an analogue of Theorem 3 for a more general model, because quantum mechanics allows us 
to consider more general quantum states and measurements. Namely, Alice can encode her message into 
a mixed state instead of a pure state and Bob can use a POVM measurement instead of an orthogonal 
measurement to recover information. A mixed state is just a probability distribution over pure states, 
so it does not provide a more general encoding model. In contrast, a POVM measurement provides a 
more general decoding model. In fact, there is another reason to extend the model. 

Example. It is not possible to construct a pure QRAC (as defined in Sect. 3.2) that simulates the 
following pure classical 2 1— > 1 RAC: 

• encoding: encode the first bit, 
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• decoding: if the first bit is requested, output the received bit; if the second one is requested — 
output no matter what is received. 



To recover the first bit with certainty, Ahce and Bob have to agree on two antipodal points on the 
Bloch sphere, where the information is encoded. Unfortunately the second bit will cause a problem — it 
is not possible to choose an orthogonal measurement of a qubit in an unknown state, so that the result 
is always the same. 

This example suggests that the model of pure quantum encoding-decoding strategies introduced in 
Sect. 3.2 should be extended in one way or the other. It is obvious that a constant decoding function (0 
or 1) can be implemented using a single-outcome POVM measurement. However, it turns out that in 
the qubit case a two-outcome POVM measurement can be replaced by a probability distribution over 
orthogonal measurements and constant decoding functions (see appendix B). This means that both 
extensions are equivalent. For simplicity we choose to extend the model by allowing constant decoding 
functions, thus Bob can either perform an orthogonal measurement or use a constant decoding function. 
The goal of this section is to show that constant decoding functions do not give any advantage. 

Definition. An enhanced orthogonal measurement is either an orthogonal measurement or one that 
always gives the same answer. 

Definition. An enhanced pure quantum n I encoding-decoding strategy is given by an ordered tuple 
{E, Ml, . . . , AIn) consisting of encoding function E : {0, 1}" i— > and n decoding functions Mi that 
are enhanced orthogonal measurements. 

Definition. An enhanced quantum encoding-decoding strategy with SR is a probability distribution over 
enhanced pure quantum strategies. 

Now it is straightforward to construct a pure quantum RAG for the previous example. In fact, now 
any classical RAG (either pure, mixed or with SR) can be simulated by the corresponding type of a 
quantum RAG. 

There is no need to further extend the model of enhanced QRAGs with SR by adding other types 
of classical randomness. For example, a probabilistic combination of POVMs does not provide a more 
general measurement, because it can be simulated by a probabilistic combination of enhanced orthogonal 
measurements. The same holds for probabilistic post-processing of the measurement results (which 
can be simulated by a probabilistic combination of enhanced orthogonal measurements as shown in 
Appendix B). Therefore enhanced QRAGs with SR constitute the most general model when any kind 
of classical randomness is allowed. 

One might suspect that the upper bound obtained in Theorem 3 does not hold for this model, but 
this is not the case. 

Theorem 4. For any n 1 enhanced QRAC with SR, p < ^ + 

Proof. According to Yao's principle and Theorem 1, we can consider the average success probability of 
pure enhanced QRAGs instead. It suffices to rule out the constant decoding functions. More precisely, 
we have to show that QRAGs having a constant decoding function for some bit give a smaller upper 
bound than those without it. In fact, we are proving a quantum analogue of Lemma 3 from Sect. 2.3.1. 

We will use induction on n. The case n = 1 is trivial — a pure enhanced QRAG with a constant 
decoding function has average success probability | < 1. Let us assume that for some n = fc — 1 > 1 
the constant decoding functions do not give any benefit. We now prove that the same holds for n = fc. 
Let us assume that the constant decoding function is used for the fcth bit. The average case success 
probability is 



x£{0,l}'' \*=i / 

where p{x, i) is the success probability (43) for the input (x, i) where i < fc— 1 and Sq^x^ is the probability 
that the decoding function gives a correct answer for the kth bit. The last bit can be ignored during 




(71) 
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the encoding and decoding of other bits: 



P(^)=|^ E 2X:p(-,^))+^ (72) 

xeio.i}''-^ i=i / 



2*^-1 • (/c- 1) ^ ' ' I 2/f 



Note that the bracketed expression in (73) is the success probabihty p{k — 1) of a shorter QRAC. 
Therefore 

p{k)^'^-p{k-l) + ^. (74) 
Now we can apply the inductive hypothesis: 



, , fc- 1 /I 1 \ 1 1 Vfc^ 11 

completing the proof. Thus the upper bound obtained in Theorem 3 holds for the general model as 
well. □ 

Observe again that for n = 2 and n = 3 this upper bound matches equations (48) and (53), 
respectively. This means that the known 2 i-^ 1 and 3 i-^ 1 quantum random access codes with pure 
encoding-decoding strategies (see Sects. 3.3.1 and 3.3.2, respectively) are optimal even among enhanced 
strategies with SR. For rt = 4 we get p < j- 

A similar upper bound was recently obtained by Ben-Aroya et al. [16] for rn^ m QRACs, where k 
bits must be recovered. They allow randomized strategies without shared randomness. In particular, 
they show that for any 77 > 2 In 2 there exists a constant C,, such that for n ^ k 

It might be possible to generalize our upper bound (75) to obtain something similar to (76). 



3.8 Lower bounds 

In the next two sections we will describe two constructions of n 1— > 1 QRAC with SR for all n > 1. 
These constructions provide a lower bound on the success probability. They use random and orthogonal 
measurements, respectively. In the first case it is hard to compute the exact success probability even 
for small values of n, but we will obtain an asymptotic expression. However, in the second case we do 
not know the asymptotic success probability, but can easily compute the exact success probability for 
small n. 



3.8.1 Lower bound by random measurements 

We now turn to lower bound for p. A lower bound for QRACs with shared randomness can be obtained 
by randomized encoding. Alice and Bob can use the shared random string to agree on some random 
orthogonal measurement for each bit. Each of these measurement bases can be specified by antipodal 
points on the Bloch sphere (see Sect. 3.1.1). These points can be sampled by using some sphere 
point picking method [21], near uniformly given enough shared randomness. The chosen measurements 
determine the optimal encoding scheme (see Sect. 3.4) which is known to both sides. 

The expected success probability of randomized n 1-^ 1 QRAC similarly to (60) is given by 

E{p)^l(l + - E d{v,,--- ,v^)] (77) 
2 V / 
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where according to equation (61) 



E d(vi, ■ ■ ■ ,Vn) = E I — V 

\ ae{l.-l} 

= ^ y E 



ae{i,-i}" 



(78) 
(79) 



Each a G {1,-1}" influences the direction of some vectors Vi, but the resulting set {uiVi} is still 
uniformly distributed. Therefore the expected value in equation (79) does not depend on a and we have 



E d{vi,--- ,i;„) = E 



(80) 



This expression has a very nice geometrical interpretation — it is the average distance traveled by a 
particle that performs n steps of unit length each in a random direction. This distance can be found 
by evaluating the following integral: 



Jip„=Q 



E 



'sin 9i cos ifi 
sin 9i sin ipi 
cos 9i 



Y[ sin 9.^ d9^ dcpi 



(81) 



Unfortunately it is very hard to evaluate it even numerically, since the integrand is highly oscillatory. 
An alternative approach is to directly simulate a random walk by sampling points uniformly from the 
sphere [21]. For small values of n the success probability (77) averaged over 10^ simulations is given in 
Table 1. Luckily, we have the following asymptotic result: 

Theorem 5 (Chandrasekhar [22, pp. 14], Hughes [23, pp. 91]). The probability density to arrive at 
point R after performing n ^ 1 steps of random walk is 



W{R) 



inn 



3/2 



exp 



2n 



(82) 



Theorem 6. For every n 3> 1 there exists an n i—s- 1 QRAC with expected success probability p 



Proof. Because of the spherical symmetry of the probability density in formula (82) , the average distance 
traveled after n ^ 1 steps of random walk is given by 



E 

From (80) and (77) we obtain 



/ R-W{R)-ATTR^dR = 2J—. 
Jo V 37r 



2 



which gives the desired lower bound. 



(83) 

(84) 
□ 



Formally this lower bound holds only for large n. However, if one estimates the actual value of (83) 
by random sampling one can see that the asymptotic expression (84) is indeed smaller than the actual 
value (see Table. 1). 



22 



3.8.2 Lower bound by orthogonal measurements 

According to the upper bound obtained in Sect. 3.6 the known 2 t-^ 1 and 3 i-^ 1 QRACs (see Sect. 3.3) 
are optimal. This suggests that orthogonal measurements can be used to construct good codes. Un- 
fortunately this idea cannot be directly applied when n > 3, since in M.^ there are only three mutually 
orthogonal directions. However, if we choose roughly one third of all measurements along each coordi- 
nate axis, we will get quite a lot of mutually orthogonal measurement pairs. 

Let Vi = (1,0,0), V2 — (0,1,0), V3 = (0,0,1), and Vi : = Vi. According to equation (60) in 
Sect. 3.5 the success probability of this n 1— > 1 QRAC with SR is related to the average distance (61) 
traveled by a random walk. For our choice of measurement directions Vi the random walk takes place 
in a cubic lattice and consists of roughly n/3 steps along each coordinate axis. Thus we can simplify 
the equation (61) for the average distance traveled to avoid having an exponential number of terms in 
it: 

d(«l, ■ • ■ . ^n) = ^ E E E (") (^) (fc) ^/(^^W^W^W^¥^\ (85) 

where x -\- y -\- z = n and each of x, y, z is roughly n/3. The corresponding success probability can be 
obtained by plugging this expression in equation (60). 

This lower bound is better than the one obtained in the previous section using random measurements 
and it also requires less shared randomness. The difference of both lower bounds is shown in Fig. 10. 
The periodic pattern of length 6 in this picture can be explained as follows. When n is a multiple of 
3, the same number of steps of a random walk is performed along each coordinate axis (this explains 
the factor 3). To explain the factor 2, let us consider a random walk on a line, i.e., one of the three 
coordinate axis. The distinction between odd an even number of steps of such a walk is that the 
probability distribution after an even number of steps is peaked at the origin, but this peak has no 
contribution whatsoever to the average distance traveled. This intuition suggests that it should be 
especially hard to beat this lower bound when n is of the form 6/c + 3. 



4 Constructions of QRACs with SR 

It is plausible that one can do better than the lower bound obtained above, which used random mea- 
surements. In this section we will consider several constructions of quantum random access codes with 
shared randomness for some particular values of n. First, in Sect. 4.1 we will describe numerically 
obtained QRACs. Then, in Sect. 4.2 we will construct new QRACs with high degree of symmetry. In 
Sect. 4.3 we will compare both kinds of codes and draw some conclusions. 



Random measurements 



Asymptotic 



0.825735 
0.765962 
0.730329 
0.706013 
0.688063 
0.674113 
0.662868 
0.653553 



Sampling 



0.8333 
0.7708 
0.7333 
0.7082 
0.6897 
0.6754 
0.6638 
0.6544 



Orthogonal measurements 



Numerical 



0.853553 
0.788675 
0.741481 
0.711803 
0.686973 
0.677458 
0.666270 
0.656893 



Exact 



+ 



— r~ 

2-^/2 

2V3 
1 + 73 
8^2 
2+V5 
20 

1 + 73+n/6 
16^3 

15+6N/5+2N/l3+yi7 
224 

12+9\/3+6\/5+6n/7+v^ 

256V^ 
1073+9^11+3719 
384 



Table 1: Comparison of n i-^- 1 QRACs with SR that use random and orthogonal measurements, 
respectively. For the first code we give the success probability according to the asymptotic expression 
(84) and a numerical value obtained by 10® random samples. For the second code we give both the 
numerical and the exact value of the success probability according to equation (85). 
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Figure 10: The difference of both lower bounds for QRACs with SR. Black squares correspond to the 
bound obtained using measurements along coordinate axes and the horizontal line corresponds to the 
asymptotic bound (84) using random measurements (see Sects. 3.8.1 and 3.8.2, respectively). The first 
bound is better, except for n = 6 (notice a periodic pattern of length 6). 



4.1 Numerical results 



n 


Section 


Probability 


2 


4.1.1 


0.853553 


3 


4.1.1 


0.788675 


4 


4.1.2 


0.741481 


5 


4.1.3 


0.713578 


6 


4.1.4 


0.694046 


7 




0.678638 


8 




0.666633 


9 


4.1.5 


0.656893 


10 




0.648200 


11 




0.641051 


12 




0.634871 



Table 2: The success probabilities of numerical n i-^ 1 QRACs. 

In this section we will discuss some particular n i-^ 1 QRACs with shared randomness for several 
small values of n. These codes were obtained using numerical optimization. The optimization must 
be performed only over all possible measurements, because in Sect. 3.4 we showed that the choice of 
measurements determines the optimal encoding in a simple way. Each measurement is specified by a 
unit vector Vi e . For n i-^ 1 QRAC there are n such vectors and one needs two angles to specify each 
of them. Without loss of generality we can assume that t^i — (0, 0, 1) due to the rotational symmetry 
of the Bloch sphere. Thus only 2(n — 1) real parameters are required to specify all Vi and therefore an 
n i-^ 1 QRAC. To find the best configuration of measurements Vi, one needs to given by 
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(63). According to (59) the success probability of the corresponding QRAC is given by 



This is not a convex optimization problem, since the feasible set (given by \\vi\\ = 1 for all 1 < i < n) 
is not convex. Note that it is not convex even if we relax equalities \\vi\\ = 1 to inequalities \\vi\\ < 1. 
We used the Mathematical & general- purpose built-in function NMaximize to solve this problem. 

Once the measurements Vi are found, one can easily obtain the Bloch vector of the qubit state 
that must be used to optimally encode the string x. We showed (see Sect. 3.4) that is a unit vector 
in direction v^, where is given by (57). For almost all QRACs that we have found using numerical 
optimization, the points r^. form a symmetric pattern on the surface of the Bloch ball. Thus we were 
able to guess the exact values of and Vi. However, as in any numerical optimization, optimality of 
the resulting codes is not guaranteed. 

In order to make the resulting codes more understandable, we depict them in three dimensions using 
the following conventions: 

• each red point encodes the string indicated, 

• each blue point defines the axis of the measurement when the indicated bit is to be output, and 

• for each measurement there is a corresponding (unlabeled) blue great circle containing states 
yielding and 1 equiprobably. 

More precisely, the blue point with label i defines the basis vector jV'o) corresponding to the outcome 
of the ith measurement (see Sect. 3.2). Note that the blue circles and blue points come in pairs — the 
vector |'(/;q) defined by the blue point is orthogonal to the corresponding circle. As a cautionary note, 
occasionally, the blue point for one measurement falls on the great circle of a different measurement 
(for example, blue points 1 and 2 in Fig. 12 lie on one another's corresponding circles). If there are too 
many red points, we omit the string labels for clarity. 

Usually the codes have some symmetry; for example, the encoding points may be the vertices of 
a polyhedron. In such cases we show the corresponding polyhedron instead of the Bloch sphere. We 
do not discuss 7 i-^ 1 and 8 i-^ 1 QRACs since the best numerical results have almost no discernible 
symmetry. We also do not discuss the numerical results for n > 10 (see Table 2 for success probabilities). 
The numerically obtained 10 i— > 1 code is symmetric and resembles 6 i-^- 1 code discussed in Sect. 4.1.4, 
but the 11 1-^- 1 and 12 i—!- 1 codes again have almost no discernible symmetry. Success probabilities of 
all numerical n^-^ 1 QRACs with SR are summarized in Table 2 and Fig. 11. 



4.1.1 The 2 ^ 1 and 3 ^ 1 QRACs with SR 

We used numerical optimization as described above to find 2 > 1 and 3 i— > 1 QRACs with shared 
randomness and obtained the optimal codes discussed in Sects. 3.3.1 and 3.3.2. 

The codes are shown in Fig. 12 and 13, respectively. In the first case the encoding points are the 
vertices of a square and the success probability is 

P=^ + ^« 0.853553. (87) 

In the second case they are the vertices of a cube. The success probability is 

p=- + « 0.788675. 
2 2\/3 



4.1.2 The 4 ^ 1 QRAC with SR 

In Sect. 3.3.3 we discussed the impossibility of a 4 i-^ 1 QRAC when Alice and Bob are not allowed to 
cooperate. However, a 4 i-^ 1 QRAC can be obtained if they have shared randomness. The particular 
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Figure 11: Success probabilities p{n) of numerical n i-^ 1 QRACs with SR from Table 2. The upper 
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respectively) . 



bound h + and the lower bound h + \ / are indicated by dashed lines (see Sects. 3.7 and 3.8.1 
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Figure 12: The 2 i-> 1 QRAC with SR." 



"For those who are using a black-and-white 
printout: this is how red and blue looks like. 
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Figure 13: The 3 1 QRAC with SR. 
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Figure 14: The 4 1 QRAC with SR. 



4 1-^ 1 QRAC with SR discussed here was found by a numerical optimization. It is a hybrid of the 
2 i-^- 1 and 3 i— > 1 codes discussed in Sects. 3.3.1 and 3.3.2, respectively. 

The measurements are performed in the bases (Mi, M2, Afa, Afa), where Mi, M2, and are the 
same as in the 3 1— > 1 case (note that the last two bits are measured in the same basis, namely M3). 
These bases are given by (44), (45), and (49), respectively. The points that correspond to an optimal 
encoding for these bases are the vertices of a regular square -^{±1, ±1, 0) in the xy plane and a cube 

^(±1,±1,±2) that is stretched in the z direction (see the Bloch sphere in Fig. 14). The Bloch vector 
for the string x — X1X2X3X4 is explicitly given by 



The encoding function can be described as follows: 

• if X3 = X4, use the usual 3 i— > 1 QRAC with an emphasis on x^ to encode the string X1X2X3, 

• if X3 7^ Xi — encode only X1X2 using the usual 2 1-^ 1 QRAC. 

In the 3 i-^- 1 scheme the probability to recover x^ must be increased by stretching the cube along the 
z axis, because X3 equals X/i and therefore it is of greater value than xi or X2- 

This 4 I— > 1 QRAC can also be seen as a combination of two 3 i— > 1 QRACs: the string 2:1X22:3 is 
encoded into the vertices of a smaller cube inscribed in a half of the Bloch ball (the vertices that lie 
within the sphere are projected to its surface). The last bit X4 indicates in which half the smaller cube 
lies (the upper and lower hemispheres correspond to 2:4 = and 1, respectively). 

The qubit state is explicitly given by E(xi, X2, 2:3, 2:4) — a\0) + f3 |1), where 




(89) 




(90) 
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The 16 values for f3 are exactly the sixteen roots of the polynomial (recall Sect. 3.1.2) 

23O4/3I6 + 3072/3^2 _,_ 1120/38 + 128/34 + 1. (91) 

If a shared random string is not available, the worst case success probability of this QRAC is |. 
However, if shared randomness is available, input randomization (as in Lemma 1) can be used and we 
will get the same success probability for all inputs, namely 

.4 + ^^^0.741481. (92) 
We do not know if this 4 1^ 1 QRAC with SR is optimal. 
4.1.3 The 5 1 QRAC with SR 

To obtain a 5 1-^ 1 QRAC, we take the bases Mi, M2, and M3, given by (44), (45), and (49), respectively, 
and also 



2 V* + 1/ ' 2 V* + 1 

The Bloch vectors V3 — (0, 0, ±1) for the basis M3 are along the z axis, but the Bloch vectors of the other 
four bases form a regular octagon in the xy plane (shown in Fig. 15): Vi = (±1,0,0), V2 = (0,±1,0), 
V4 = ±^(1, 1,0), — ±^(—1, 1,0). The Bloch vector encoding the string x = a;ia;2X3a;4X5 is 

^ /\/2(-l)-i+(-l)-^-(-l)-A 
r{x) = \/2(-l)"^ + (-1)"^ + (-1)"^ , (95) 

y^lO + s(x)4^/2 \ V2(-l)"^ / 

where s{x) £ {—1, 1} and is given by 

s{x) = ^ ^ ^ (-1)-- - ^ ^ / ^ (-1)-^ (96) 
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1 

Figure 16: The "preferable regions" of the measurement M2 (only the upper hemisphere is shown, 
the other half is symmetric). For each of the measurements the direction of the Bloch vector {ipo) is 
indicated by the corresponding number. The white triangles correspond to X2 = 0, but the gray ones 
to a;2 = 1. 



The great circles with equiprobable outcomes of the measurements partition the Bloch sphere into 
16 equal spherical triangles. There are two strings encoded into each triangle. The idea for how to 
locate the correct point for the given string x is as follows. Observe that the strings with X3 = and 
X3 — 1 are encoded into the upper and lower hemisphere, respectively (this means that for all strings 
the probability that the measurement M3 gives the correct value of X3 is greater than ^ ) . Next observe 
that half of all strings have s{x) = 1, but the other half have s{x) = — 1 (in fact, the two strings in the 
same triangle have distinct values of s). 

Let us first consider the case s{x) = 1. We call such string compatible with the measurements, 
because it can be encoded in such a way that every measurement gives the correct value of the corre- 
sponding bit with probability greater than ^ . For the ith bit of x we can define the "preferable region" 
on the Bloch sphere as the hemisphere where Mi recovers Xi with probability greater than i. The 
intersection of these five regions is one sixteenth of the Bloch sphere — the triangle where x must be 
encoded. The point with the smallest absolute value of the z coordinate in this triangle must be chosen 
(it has smaller probability than other points in the triangle to recover X3 correctly, but the probabilities 
for the other four bits are larger). 

If s{x) = — 1, the string x is incompatible with the measurements, because the intersection of the 
"preferable regions" is empty. Thus, no matter where the string is encoded, at least one bit will differ 
from the most probable outcome of the corresponding measurement. We can take this into account and 
modify the definition of the "preferable region" for the ith bit (i ^ 3). It is a union of eight triangles: 
four triangles where the most probable outcome of Mi equals Xi , and four triangles where it does not 
equal Xi (in either case the triangles with maximal probability of correct outcome of Mi must be taken) . 
For example, the "preferable regions" for X2 are shown in Fig. 16. The regions for X3 remain the same 
as in the previous case. The intersection of all five regions for the given string x is the triangle where 
the string must be encoded. The point in the triangle with the largest absolute value of the z coordinate 
must be chosen. As a result, three of the measurements will give the correct value of the corresponding 
bit of the string x with probability greater than ^. 

The corresponding qubit state is given by E{xi, X2, X3, X4, x^) = a \ 
(3 defined as follows: 



- f3\l) with coefficients a and 



(-1)^ 



/3 = 



(-l)-i+z(-l)-^ + ^(- 



-l)-^ + ^(-l) 



(97) 



10 + s{x)AV2 + 2(-l)^3 ^5 + s{x)2V2 
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Figure 17: The measurements for the 6 ^ pigyre 18: The 6 ^ 1 QRAC with SR. 

1 QRAC shown on the right. 



The coefhcients /3 are the roots of the polynomial 

1336336/?^^ + 961792/324 + 15U32(3^^ + 1600/3^ + 1. (98) 
Again, using input randomization we obtain the same success probability for any input, namely 



P^l + ^ \/2{5 + Vl7) « 0.713578. (99) 



4.1.4 The 6 ^ 1 QRAC with SR 

The Bloch vectors corresponding to the 6 measurements are as follows: 



t;i =±(0,+l,+l)/V2, 
t;2 = ±(0,-l,+l)/\/2, 
t;3-±(+l,0,+l)/V2, 



(100) 

= ±(+l,0,-l)/\/2, 



Ve = ±{-l,+l,0)/V2. 

They correspond to the 12 vertices of the cuboctahedron (or the midpoints of the 12 edges of the cube) 
and are shown in Fig. 17. The great circles orthogonal to these vectors form the projection of the edges 
of a normalized^ tetrakis hexahedron and partition the Bloch sphere into 24 parts (see Fig. 18). Each 
of these parts contains one vertex of a truncated octahedron — the dual of tetrakis hexahedron. It is 
inscribed in the Bloch sphere shown in Fig. 18. 



®The vertices of the tetrakis hexahedron are not all at the same distance from the origin (the ones forming an octahedron 
are 2/V3 times closer than those forming a cube). So the polyhedron has to be normalized to fit inside the Bloch sphere 
(the vectors pointing to the vertices have to be rescaled to have a unit norm). 
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The measurement bases corresponding to Vi can be found using (33): 



Ml = 

Mi = 
Ms = 

M6 = 



1 / V2 + \/2\ 1 / V2-\/2 ^ 

2 l^iV2- V2; ' 2 l^-iV2 +\/2^ 

1 / y^2 + 72 \ l( V2-V2^ 

2 -V2J ' 2 l^zv/2 + x/2^ 

1 / V2 + V2\ 1 / V2-V2 \ ' 

2 Iy2^r2 l_y2Wl/ 



(101) 



1 / V2-V2\ 1 / V2 + V2 ^ 

2 l^v/2+72; '2 \-V2^72j 

l/^/2^ 1/ ^/2 
2 V+V '2 V-«-l 
l/y2\ If V2 
2 U- 1/ ' 2 Ui + 1 



Note that and A/g are the same as (93) and (94) for the 5 1-^ 1 QRAC described in the previous 
section. Another way to describe these 6 bases is to consider the (3 coefficients for the 12 vectors that 
form them. It turns out that these coefficients are exactly the roots of the polynomial 



256/3^2 - 128/3^ - 44/3^^ + 1. 



(102) 



Let us consider how to determine the point where a given string should be encoded. According to 
(57) we have to find the sum of the vectors Vi defined in (100), each taken with either a plus or a minus 
sign. These vectors correspond to six pairs of opposite edges of a cube and the signs determine which 
edge from each pair we are taking (see Fig. 17). There are only three distinct ways of doing this (see 
Fig. 19). Regardless of which way it is, for each of the chosen edges there is exactly one other that 
shares a common face and is parallel to it. Thus we can partition the chosen edges into three pairs (in 
Fig. 19 such pairs are joined with a thick blue line). The sum of the vectors Vi for edges in a pair is 
always parallel to one of the axes and its direction is indicated with an arrow in Fig. 19. From these 
arrows one can see where the encoding point should lie. 

Now let us classify all 2^ = 64 strings of length 6 into 3 types according to the location of the encoding 
point on the Bloch sphere. Each type of string is encoded into a vertex of a specific polyhedron (see 
Fig. 20). These polyhedra are the cube, the truncated octahedron, and the octahedron and the number 
of strings of each type are 16, 24, and 24, respectively. Let us consider them case by case: 



Truncated 
octahedron 


Octahedron 


* * 1110 


* * 1101 


* * 0001 


* * 0010 


10 * *11 


01 * *11 


01 * *00 


10 * *00 


1110 * * 


1101 * * 


0001 * * 


0010 * * 



Table 3: Patterns of strings corresponding to the vertices of truncated octahedron and octahedron ("*" 
stands for any value). 



• The cube has 8 vertices: 



V3 



(±1,±1,±1) 



(103) 
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Figure 19: Three distinct ways of choosing one edge from each pair of opposite edges of a cube. The 
chosen edges are marked with blue points. Points lying on opposite edges of the same face are connected 
and the direction of the sum of the corresponding vectors is indicated with an arrow. The corresponding 
encoding point is shown in red. The red points obtained from all possible choices of the same kind are 
the vertices of a cube, a truncated octahedron, and an octahedron, respectively (see Fig. 20). 



and there are 2 strings encoded into each vertex. These 16 strings are exactly those xiX2 ■ . - Xe S 
{0, 1}^ that satisfy 

-X2\ + \x3- X4\ + \X5 - xel e {0, 3} . (104) 
This condition ensures that the three arrows in Fig. 19 are orthogonal. 

• The truncated octahedron has 24 vertices. Their coordinates are obtained by all permutations of 
the components of 

-^(0,±1,±2). (105) 

There is just 1 string encoded into each vertex. In this case there will be two pairs of chosen 
edges that belong to the same face (note the "cross" in the Fig. 19 formed by pairs whose arrows 
are pointing outwards of the page). The third pair (with the arrow pointing up) can be rotated 
around this face to any of the four possible positions. This corresponds to fixing four bits of the 
string and choosing the remaining two bits in an arbitrary way. Since the "cross" can be on any 
of the six faces of the cube, one can easily describe all 24 strings of this type (they are listed in 
the first column of Table 3). 

• The octahedron has 6 vertices: 

(±1,0,0) U (0,±1,0) U (0,0,±1) (106) 

and there are 4 strings encoded into each vertex. In this case two arrows in Fig. 19 are pointing to 
opposite directions (up and down). If we fix these arrows, we can rotate the third one (pointing 
outwards) in any of four directions. Hence we can describe all 24 strings of this type in a similar 
way (see the second column of Table 3). 

The coefficients f3 of the encoding states are the 64 roots of the polynomial 

P^iP - 1)4(4/3* - 1)4(36/3^ + 24/3* + 1)^ 
(25(3^ - 15/3* + 1) (400/3^ - 360/3* + 1) (400/3^ + 56/3* + 25). (107) 
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Figure 20: Three polyhedra (cube, truncated octahedron, and octahedron) corresponding to three differ- 
ent types of strings for 6 i-^- 1 QRAC with SR. The red points in Fig. 18 are obtained by superimposing 
these three polyhedra. 



The obtained success probabihty using input randomization is 

1 2 + x/3 + %/T5 , , 
P = 2 + « 0.694046. (108) 

4.1.5 The 9 ^ 1 QRAC with SR 

This QRAC is a combination of three 3 i-^ 1 QRACs described in Sect. 3.3.2. It has three measurements 
along each axis: 

vi = Vi = Vj = ±(1,0,0), 

T^2 = t^5 = ^^8 = ±(0,1,0), (109) 

vj,^ Vq = Vq = ±(0,0, 1). 

The measurement bases Mi, M2, and M3 corresponding to the Bloch vectors V2, and are given 
by (44), (45), and (49), respectively. 

The encoding points can be characterized as a 4 x 4 x 4 cubic lattice formed by vectors (57) projected 
on the surface of the Bloch ball. Note that this lattice consists of vertices of 8 equal cubes each lying 
in a different octant. Then the 7 points inside of each spherical triangle in Fig. 21 are the projection of 
the vertices of the corresponding cube.^ 

All 2^ = 512 strings can be classified into 3 types. First consider a string 010203 € {0, 1}"^ and define 

|ai - 02! + |a2 - 03I + I03 - Oil 
5(01,02,03) = ^ . (110) 

Notice that 5(01,02,03) <E {0, 1}. Now for x = X1X2 . . . xg € {0,1}^ define 

t{x) = s(a;i,a;4,a;7) + s(x2,X5,a;8) + 5(0:3, Xg, Xg). (Ill) 
Then the type of the string x can be determined as follows: 

0, 3 cube, 

t{x) — {1 truncated cube, (112) 

2 small rhombicuboctahedron. 

These types are named after polyhedra, since each type of string is encoded into the vertices of the 
corresponding polyhedron (see Fig. 22): 



We get 7 points instead of 8 since the projections of two diagonally opposite vertices coincide. 
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Figure 21: The 9 1 QRAC with SR. 



• The cube has 8 vertices and there are 28 strings encoded into each vertex. These vertices are: 

-^(±1,±1,±1). (113) 

• The deformed® truncated cube has 24 vertices and there are 3 strings encoded into each vertex. 
These vertices are: 

^(±1, ±3, ±3) U ^(±3, ±1, ±3) U ^(±3, ±3, ±1). (114) 

• The deformed^ small rhombicuboctahedron also has 24 vertices and there are 9 strings encoded 
into each vertex. These vertices are: 

^(±3,±1,±1) U -L(±l,±3,±l) U ^(±1,±1,±3). (115) 

Vn Vn Vn 

The coefficients /3 for the corresponding qubit states a\0) + (3\1) are the roots of the following 
polynomial: 

(36/3** + 24^^ + 1)28(1444/3* + 760/3"* + 81)^(484/3** + 440/3"* + 1)^ 
(52128400/3*^ - 21509824/3*^ + 26780424/3® - 372400/3* + 15625)^ 
(5856400/3*^ - 1788864/3*^ + 1232264/3® - 92400/3* + 15625)^ (116) 

Using input randomization we get success probability 



1 , 192+10^/3 + 9^/n + 3^A9 ^^^^^^^^ 

2 384 0.656893. (117) 



*The edges of the truncated cube are of the same length. In our case the eges forming triangles are \/2 times longer 
than the other edges. 

^The edges of the small rhombicuboctahedron are also of the same length, but in our case the edges forming triangles 
again are \/2 times longer. 
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Figure 22: Three polyhedra (cube, small rhombicuboctahedron, and truncated cube) corresponding to 
three different types of strings for 9 i-^ 1 QRAC with SR. The red points in Fig. 21 are obtained by 
superimposing these three polyhedra. 



4.2 Symmetric constructions 

In Sect. 4.1 we have discussed in great detail n i-^ 1 quantum random access codes with shared ran- 
domness for some particular values of n. Since these codes were obtained using numerical optimization, 
there are still some questions left open. Most importantly, are the codes for n > 4 discussed in Sect. 4.1 
optimal? If this is the case, do these codes (see Figs. 12, 13, 14, 15, 18, and 21) have anything in 
common that makes them so good? 

The purpose of this section is to shed some light on these two questions. We will explore the 
possibility that symmetry is the property that makes QRACs with SR good. In Sect. 4.2.1 we will 
explore what symmetries the codes found by numerical optimization have and what other symmetries 
are possible. In several subsequent sections we will use these symmetries to construct new codes and 
compare them with the numerical ones (the success probabilities of the obtained codes are summarized 
in Table 4). In Sect. 4.3 we will conclude that symmetric codes are not necessarily optimal and speculate 
about what else could potentially be used to construct good QRACs. 



n 


Section 


Probability 


4 


4.2.2 


0.733253 


6 


4.2.3 


0.694042 


9 


4.2.4 


0.656393 


15 


4.2.5 


0.620183 



Table 4: The success probabilities of symmetric n i-^ 1 QRACs with SR. See Table 7 for the comparison 
with numerically obtained codes. 



4.2.1 Symmetric great circle arrangements 

If we want to construct a QRAC with SR that has some sort of symmetry, we have to choose the 
directions of measurements in a symmetric way. In other words, we have to symmetrically arrange the 
great circles that are orthogonal to the measurement directions. 

In this section we will discuss two ways that great circles can be arranged on a sphere in a symmetric 
way. These arrangements come from quasiregular polyhedra and triangular symmetry groups, respec- 
tively. The first kind of arrangement is not directly observed in numerically obtained examples, despite 
its high symmetry. However, the second one is observed in almost all numerically obtained codes. Since 
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Figure 23: Quasiregular polyhedra: cuboctahedron and icosidodecahedron. 



our approach is empirical, we will not justify when an arrangement is "symmetric enough" to be of 
interest. We will use the term symmetric codes to refer to the codes constructed below. This is just 
to distinguish them from numerically obtained codes in Sect. 4.1, not because they satisfy some formal 
criterion of "being symmetric" . 

Quasiregular polyhedra 

A (convex) quasiregular polyhedron is the intersection of a Platonic solid with its dual. There are 
only three possibilities: 



The tetrahedron is self-dual thus the octahedron, which is the intersection of two tetrahedrons, has 
slightly different properties than the other two polyhedra (e.g., its all faces are equal). For this reason 
octahedron may be considered as a degenerate quasiregular polyhedron or not be considered quasiregular 
at all since it is Platonic. Thus there are only two (non-degenerate) convex quasiregular polyhedra (see 



These polyhedra have several nice properties. For example, all their edges are equivalent and there 
are exactly two types of faces (both regular polygons), each completely surrounded by the faces of 
the other type. The most relevant property for us is that their edges form great circles. Since the 
arrangements of great circles formed by the edges of cuboctahedron and icosidodecahedron do not 
appear in the numerical codes, we will use them in Sects. 4.2.2 and 4.2.3 to construct new (symmetric) 
4 I— !■ 1 and 6 i-^ 1 QRACs with SR, respectively. 

Triangular symmetry groups 

Consider a spherical triangle — it is enclosed by three planes that pass through its edges and the 
center of the sphere. Let us imagine that these planes are mirrors that reflect our triangle. These three 
reflections generate a reflection group [24, 25]. For some specific choices of the triangle this group is 
finite and the images of the triangle under different group operations do not overlap. Hence they form 
a tiling of the sphere. This tiling can also be seen as several (most likely more than three) great circles 
cutting the sphere into equal triangles. 

We can choose any of the triangles in the tiling and repeatedly reflect it along its edges so that it 
moves around one of its vertices. This means that the angles of the corners that meet at any vertex of 

^"Several possible criteria are: (a) any great circle can be mapped to any other by a rotation from the symmetry group 
of the arrangement, (b) the sphere is cut into pieces that are regular polygons, (c) the sphere is cut into pieces of the 
same form. However, not all examples we will give satisfy these three conditions. In fact, each condition is violated by at 
least one of the examples we will consider. 



cuboctahedron — cube n octahedron, 
icosidodecahedron — icosahedron n dodecahedron. 



octahedron = tetrahedron H tetrahedron 



(118) 
(119) 
(120) 



Fig. 23). 
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Figure 24: Triangular symmetry groups. First row: (2,2,2), (2,2,3), (2,2,4). Second row: (2,3,3), 
(2,3,4), (2,3,5). 

the tiling must be equal. Moreover, we do not want the triangle to intersect with any of the mirrors, 
so only an even number of triangles can meet at a vertex. 

Hence the angles of the spherical triangle must be for some integers p,q,r > 2. The sum 

of the angles of a spherical triangle is at least tt, so the numbers p, q, r must satisfy 

- + - + ->!. (121) 
p q r 

lip <q<r, the only solutions are: (2,2,fc) for any fc > 2, (2,3,3), (2,3,4), and (2,3,5). The tilings 
corresponding to these solutions are shown in Fig. 24. The symmetry group of such tiling is called 
triangular symmetry group [25, pp. 158] and is denoted by {p,q,r). 

We can observe these tilings in almost all numerically obtained QRACs discussed in Sect. 4.1. They 
are formed when the great circles corresponding to measurements partition the Bloch sphere into equal 
triangles. All such cases are summarized in Table 5. Tilings appearing in 2 i-^ 1 and 4 i-^ 1 QRACs 
that are not mentioned in the table can be seen as degenerate cases. 



n 


{p, r) 


Polyhedron 


Section and fig 


;ure 


3 


(2,2,2) 


octahedron 


Sect. 4.1.1, Fig 


5. 13 


5 


(2,2,4) 


normalized octagonal dipyramid 


Sect. 4.1.3, Fig 


5. 15 


6 


(2,3,3) 


normalized tetrakis hexahedron 


Sect. 4.1.4, Fig 


5. 18 


9 


(2,2,2) 


octahedron 


Sect. 4.1.5, Fig 


21 



Table 5: Triangular symmetry groups of numerical n i— > 1 QRACs. 



The tilings corresponding to triangular symmetry groups (2,3,4) and (2,3,5) do not appear in nu- 
merically obtained codes. Thus we will use them to construct new (symmetric) 9 i— > 1 and 15 1 
QRACs with SR in Sects. 4.2.4 and 4.2.5, respectively. To each tiling one can associate a corresponding 
polyhedron with equal triangular faces. The polyhedra corresponding to tilings (2, 3, 4) and (2, 3, 5) are 

^^Fore example, if we project the edges of an icosahedron on the sphere, wo obtain arcs that form a tiUng with five 
triangles meeting at each vertex. We cannot use these arcs as mirrors, since they do not form great circles (we cannot 
extend any of them to a great circle, without intersecting other triangles). 
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0011 

Figure 25: Symmetric 4 i-^ 1 QRAC with SR. 



called the normalized^^ disdyakis dodecahedron and the normalized disdyakis triacontahedron, respec- 
tively. 

Polyhedra arising from both types of symmetric great circle arrangements (quasiregular polyhedra 
and triangular symmetry groups) are summarized in Table 6. The great circle arrangements corre- 
sponding to the four marked polyhedra do not appear in numerically obtained codes, so we will use 
them to construct new (symmetric) QRACs with SR. 



n 


Faces 


(P, r) 


Polyhedron 


3 


8 


8 


(2,2,2) 


octahedron 


4 


14 


14 


QR 


cuboctahedron / 


6 


32 


32 


QR 


icosidodecahedron / 


6 


24 


32 


(2,3,3) 


normalized tetrakis hexahedron 


9 


48 


74 


(2,3,4) 


normalized disdyakis dodecahedron / 


15 


120 


212 


(2,3,5) 


normalized disdyakis triacontahedron / 



Table 6: Polyhedra whose edges form great circles. The first column indicates the number of great 
circles. The next two indicate, respectively, the number of faces of the polyhedron and the maximal 
number of pieces achievable by cutting the sphere with n great circles (see Sect. 3.3.3). The fourth 
column indicates the triangular symmetry group (QR means quasiregular). The name of the polyhedron 
is given in the last column. Four marked polyhedra will be used in subsequent sections to construct 
symmetric QRACs with SR. 



4.2.2 Symmetric 4 ^ 1 QRAC with SR 

Recall that in Sect. 3.3.3 we proved that four planes passing through the center of the Bloch sphere 
partition its surface into at most 14 parts. The most symmetric way to obtain 14 parts is to use the 
four planes parallel to the four faces of a regular tetrahedron. The measurements are along the four 
directions given by the vertices (see Fig. 26). 

The simplest way to construct a regular tetrahedron is to choose four specific vertices of the cube, 
i.e., from the set ^(±1, ±1, ±1). For example, we could choose the ones with an odd number of positive 
coordinates. They provide us with the following pairs of antipodal Bloch vectors as the measurement 

Normalized means that all vectors pointing from the origin to the vertices of the polyhedron are rescaled to have unit 
norm. 



38 



Figure 26: A regular tetrahedron and four great circles parallel to its faces. The circles are determined 
by the measurements in the direction of the vertices of the tetrahedron. The numbers at the vertices 
indicate the Bloch vectors of basis states 1-00) of the measurements for the 4 i— > 1 QRAC shown in 
Fig. 25. 



bases: 

vs 

The qubit states corresponding to these Bloch vectors are as follows: 

Ml = Af(+1,+1), 
M2 = Af(+1,-1), 

(123) 

M3 = M(-1,+1), 
M4 = M(-1,-1), 

where 

The great circles determined by these measurements partition the Bloch ball into 14 parts. In 
fact, the grid formed by these circles is a projection of the edges of a cuboctahedron (see the part on 
quasireglar polyhedra in Sect. 4.2.1) on the surface of the Bloch ball (see Figs. 25 and 26). 

In each of the 14 parts of the Bloch sphere a definite string can be encoded so that each bit can be 
recovered with a probability greater than i. Strange as it may seem, the remaining 2 strings (cc = 0000 
and X = 1111) can be encoded anywhere without affecting the success probability of this QRAC. This 
is not a surprise if we recall from Sect. 3.4 that the optimal encoding of the string a; is a unit vector 
in the direction of given by equation (57). In our case the Bloch vectors of the measurement bases 
point to the vertices of a regular tetrahedron centered at the origin. They clearly sum to zero, so 
'I'oooo = ■Willi = 0- Thus the scalar product ■ v,^. in (56) is also zero and the success probability does 
not depend on the vectors roooo ^^id rmi. Therefore, we will ignore these two strings in the following 
discussion. 

The other 14 strings are encoded into the vertices of a normalized tetrakis hexahedron (the convex hull 
of the cube and octahedron). The string x = xiX2X3Xi is encoded into the Bloch vector r(x) = r^(a;), 
where 

w = xi® X2® X3® X4 e {0,1} (125) 



±(+l,-l,-l)/V3, 
±(-l,+l,-l)/\/3, 
±(-l,-l,+l)/V3, 
±(+l,+l,+l)/\/3. 



(122) 
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Figure 27: The relationship between the strings encoded into the spherical square and the adjacent 
spherical triangles according to the 4 i— > 1 QRAC shown in Fig. 25. 



is the parity of the input. In the case w = the encoding points are the vertices (±1, 0, 0) U (0, ±1, 0) U 
(0,0, ±1) of an octahedron: 



Note that the Bloch vectors ri{x) are the vertices of the same cube as the Bloch vectors of the 3 i~> 1 
QRAC discussed in Sect. 3.3.2. 

One can observe the following properties of this encoding. The surface of the Bloch ball is partitioned 
into 6 spherical squares and 8 spherical triangles. Strings with w = and w = 1 are encoded into squares 
and triangles, respectively. If w = 1 {x = 1000 or a; = 0111 and their permutations), the string has one 
bit that differs from the other three. Such a string is encoded into the basis state of the corresponding 
measurement so that this bit can be recovered with certainty. If w — 0, the string is encoded into a 
square and has the following property: each of its bits takes the value that occurs more frequently at 
the same position in the strings of the four neighboring triangles (see Fig. 27 as an example) . 

The corresponding encoding function is E{x) = Ofu, |0) + /3u, |1) with coefficients ao, /3o ^^nd ai, f3i 
explicitly given by 



\1 - - X4\ 

But for ui = 1 we get the vertices (±1, ±1, ±1)/a/3 of a cube: 




(126) 




(127) 




X3X4 + (-1) 




1 - |a;3 - X4\ 
2 

^ 1 - |xi - X4:\ + i{l - \X2 - Xi\) 



V2 



(128) 



and 




(129) 




(130) 



I3{I3 - l)(4/34 - l)(36/3^ + 24/?^ + 1). 



(131) 
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Figure 28: Symmetric 6 i-^ 1 QRAC with SR. 



Using input randomization we get the same success probabihty for any input: 

1 2 + V3 , , 

p = - + w 0.733253. (132) 

It is surprising that despite higher symmetry (compare Figs. 14 and 25) this QRAC has a lower success 
probabihty than the 4 1 QRAC discussed in Sect. 4.1.2. 

4.2.3 Symmetric 6 i-> 1 QRAC with SR 

According to the discussion in Sect. 3.3.3, six great circles can cut the sphere into at most 32 parts. 
It turns out that there is a very symmetric arrangement that achieves this maximum. Observe that 
the dodecahedron has 12 faces and diametrically opposite ones are parallel. For each pair of parallel 
faces we can draw a plane through the origin parallel to both faces. These six planes intersect the 
sphere in six great circles that define our measurements. They are the projections of the edges of the 
icosidodecahedron (see Fig. 23), which is one of the quasiregular polyhedra discussed in Sect 4.2.1. 

There is another way to describe these measurements. Notice that the icosahedron (the dual of the 
dodecahedron) has 12 vertices that consist of six antipodal pairs. Our measurements are along the six 
directions defined by these pairs. The coordinates of the vertices of the icosahedron are as follows: 

^=i=(o,±T,±i) u -=i=(±i,o,±T) u -=i=(±T,±i,o), (133) 

where r = is the golden ratio (the positive root of a;^ = a; + 1). 

Each of the 64 strings is encoded either in a vertex of an icosahedron or dodecahedron. They have 
12 and 20 vertices, respectively, so there are two strings encoded in each vertex. The union of the 
icosahedron and the dodecahedron is called the pentakis dodecahedron (see the polyhedron in Fig. 28). 

The success probability of this code is 

P = ^ + || + ^\/75 + 30V5« 0.694042. (134) 

4.2.4 Symmetric 9 1 QRAC with SR 

This code is based on the triangular tiling of the sphere whose symmetry group is (2,3,4). The great 
circles corresponding to measurements coincide with the projection of the edges of the normalized 
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disdyakis dodecahedron. We can think of this QRAC as the union of 3 i— > 1 and 6 i— > 1 codes. The first 
three measurements are along the coordinate axis as in the 3 1 QRAC discussed in Sect. 3.3.2. The 
remaining six measurements are exactly the same as for the 6 i— > 1 code discussed in Sect. 4.1.4 (see 
Figs. 17 and 18), i.e., they are along the six antipodal pairs of 12 vertices of the cuboctahedron shown 
in Fig. 23. Note that a great circle of the first kind cannot be transformed to a great circle of the second 
kind via an operation from the symmetry group of the code.^'^ 

The resulting QRAC is shown in Fig. 29 and its success probability is 

p « 0.656393. (135) 
4.2.5 Symmetric 15 1 QRAC with SR 

The triangular symmetry group of this code is (2, 3, 5) and the great circles coincide with the projection 
of the edges of the normalized disdyakis triacontahedron. To understand what the measurements are in 
this case, note that the icosidodecahedron (see Fig. 23) has 30 vertices. Their coordinates are: 

(±1,0, 0) U (0, ±1, 0) U (0, 0, ±1), (136) 
^(±1, ±T, ±t2) U ^(±t2, ±1, ±r) U ^(±r, ±T^ ±1). (137) 

The measurement directions are given by 15 antipodal pairs of these vertices. 
The obtained QRAC is shown in Fig. 30. Its success probability is 

p« 0.620183. (138) 

4.3 Discussion 

In this section we will compare and analyze the numerical and symmetric QRACs with SR described 
in Sects. 4.1 and 4.2, respectively. Hopefully these observations can be used to find new n i-^- 1 QRACs 
with SR or to generalize the existing ones (see Sect. 5.3 for possible generalizations). 

The success probabilities of numerical and symmetric QRACs with SR are given in Tables 2 and 4, 
respectively (see Table 7 for the comparison) . We see that none of the symmetric codes discussed in 
Sect. 4.2 is optimal. However, the success probabilities of numerical and symmetric codes do not differ 
much. Moreover, recall that there are two more symmetric codes (3 i— > 1 and 6 i— > 1) that coincide with 

^^For the other three symmetric codes we can transform any circle to any other in this way, i.e., the symmetry group 
acts transitively on the circles. 



42 



the numerically obtained ones (see Table 6) . Concerning these two codes we can reach more optimistic 
conclusions: the 3 i-^ 1 QRAC is optimal (see Sect. 3.6) and possibly the 6 1 QRAC (see Sect. 4.1.4) 
is as well, since we did not manage to improve it in Sect. 4.2.3. 



n 


Section 


Probability 




4.1.2 


0.741481 


4 


4.2.2 


> 0.733253 


6 


4.1.4 


0.694046 


4.2.3 


> 0.694042 




4.1.5 


0.656893 


9 


4.2.4 


> 0.656393 


15 




0.620355 


4.2.5 


> 0.620183 



Table 7: Comparison of the success probabilities of n i-^ 1 QRACs with SR. For each n the first 
probability corresponds to a numerical code, but the second one to a symmetric code. For n = 15 we 
do not have numerical results, so we just use five measurements along each coordinate. In fact, the 
numerical 4 i-^ 1 and 9 i— > 1 QRACs also use measurements only along coordinate axis. The 6 i-^ 1 
QRAC with two measurements along each coordinate axis has success probability 0.686973. 

We just saw that symmetric QRACs are not necessarily optimal. One could ask if there are other 
heuristic methods that potentially could be used to construct good QRACs with SR. We will give a 
few speculations in the remainder of this section. In particular, we will discuss some special kinds of 
measurements that could be useful. To make the discussion more general, we will not restrict ourselves 
to the case of a single qubit. 

Definition. Two orthonormal bases Bi and B2 of C' are called mutually unbiased bases (MUBs) if 
I (■(/'1 1-02)1^ = 2 ^'-'^ IV-"!) € Bi and \4>2) G B2- The maximal number of pairwise mutually unbiased 
bases in is d + 1. [26] 

When d — 2, equation (37) implies that Bloch vectors corresponding to basis vectors of different 
mutually unbiased bases are orthogonal^''. There are three such bases in and their Bloch vectors 
correspond to the vertices of an octahedron. For example, the bases Mi, M2, and M3 defined in 
Sects. 3.3.1 and 3.3.2 are MUBs (they correspond to measuring along x, y, and z axis). 

Note that the measurements for numerical 2i-^l,3n-l,4f-^l, and 9 1— s- 1 QRACs are performed 
entirely using MUBs and three out of five measurement bases for numerical 5 1 QRAC are also 
MUBs. 

There is another very special measurement that appears in our QRACs. 

Definition. A set of d^ unit vectors G is called symmetric, informationally complete POVM 
(SIC-POVM) if 1(0. IV,) I' = for any [27] 

For d = 2 there are four such quantum states. Again, from equation (37) we see that the inner 
product between any two Bloch vectors corresponding to these states is — ^. Such equiangular Bloch 
vectors are exactly the vertices of a tetrahedron, e.g., Vx, V2, fs, V4 defined in (122). They were used 
in Sect. 4.2.2 to construct a symmetric 4 1-^ 1 QRAC. 

Let us compare numerical and symmetric 4 1-^ 1 QRACs from Sects. 4.1.2 and 4.2.2, respectively. 
The first one is based on MUBs and is not very symmetric. Moreover, it looks like we are wasting 
one out of four bits, since two measurements are along the same direction. However, all measurement 
directions in the Bloch sphere are mutually orthogonal, except the ones that coincide. The second 4 i-^- 1 
code is based on a SIC-POVM and is very symmetric. However, it appears that in this case we are 
wasting two out of 16 strings, since the way we encode them does not influence the success probability. 

^^The notion of the Bloch vector can be generalized for d > 2 (see [28]). Then a similar duality holds as well (see 
equation (142) in Sect. 5.3): mutually unbiased quantum states correspond to orthogonal Bloch vectors, but orthogonal 
quantum states correspond to "mutually unbiased" Bloch vectors, i.e., equiangular vectors pointing to the vertices of a 
regular simplex. 
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Now, if we compare the success probabilities of both 4 i-^ 1 codes (see Table 7), we see that the first 
one is clearly better. Hence we conclude that 

orthogonality of the measurement Bloch vectors 
seems to be more important than symmetry. 

One can come to a similar conclusion when comparing 9 i— > 1 and 15 > 1 codes. Thus it looks like 
using roughly n/3 measurements along each coordinate axis is quite a good heuristic for constructing 
n^l QRAC with SR (see Sect. 3.8.2). 



5 Conclusion 
5.1 Summary 

We study the worst case success probability of random access codes with shared randomness. Yao's 
principle (see equation (5) in Sect. 2.2) and input randomization (see Theorem 1) is applied to consider 
the average case success probability instead (this works in both classical and quantum cases). 

In Sect. 2.3.2 we construct an optimal classical n 1 RAG with SR as follows (see Theorem 2): 
Alice XORs the input string with n random bits she shares with Bob, computes the majority and sends 
it to Bob; if the ith bit is requested. Bob outputs the ith bit of the shared random string XORed with 
the received bit. The asymptotic success probability of this code is given by equation (27) in Sect. 2.3.2: 

pH~^ + ^=- (139) 

The worst case success probability of an optimal quantum RAG with SR satisfies the following 
inequalities: 

^ + /^^^(-)4 + 27^- 

These upper and lower bounds are obtained in Sects. 3.6 and 3.8.1, respectively. 

The success probabilities of classical and quantum RAGs with SR are compared in Fig. 31. 



5.2 Open problems for n 1 QRACs 

Lower bound by orthogonal measurements. The known 2 s- 1 and 3 i-^ 1 QRAGs (see Sect. 3.3) and 
our numerical 4 i— > 1 and 9 i— > 1 QRAGs with SR (see Sects. 4.1.2 and 4.1.5) suggest that MUBs can be 
used to obtain good QRAGs (see Sect. 4.3). Indeed, m-^ 1 QRAG with orthogonal measurements (see 
Sect. 3.8.2) is better than the one with random measurements (see Sect. 3.8.1). However, we were not 
able to obtain an asymptotic expression for its success probability. This is equivalent to obtaining an 
asymptotic expression for (85), i.e., the average distance traveled by a random walk with roughly n/3 
steps along each coordinate axis. 

In Fig. 32 we show how close both lower bounds and the success probabilities of numerical QRAGs 
are relative to the upper bound from Sect. 3.6. Assume that Alice and Bob are given a point in the 
light gray region in Fig. 32 and asked to construct a QRAG with SR whose success probability is at 
least as good. Then they can use measurements along coordinate axis as in Sect. 3.8.2. If the point is 
in the dark gray region, they can use one of the numerical codes from Sect. 4.1. However, if it is in the 
white region, they have to solve the next open problem. 

Optimality of num,erical codes. Prove the optimality of any of the numerically obtained n i-^- 1 
QRAGs with SR for n > 4 discussed in Sect. 4.1. Arc the optimal constructions unique (up to isomor- 
phism)? 

Prove the "Homer conjecture" that quantum RAGs with SR are at least as good as their classical 
counterparts in the sense discussed at the end of Sect. 3.4. 
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Figure 31: Comparison of success probabilities of classical and quantum RACs. Black dots correspond 
to optimal classical RACs and the dotted line shows the asymptotic behavior. Circles correspond to 
numerical QRACs and dashed lines to quantum upper and lower bounds, respectively. 



5.3 Possible generalizations 

There are several ways that random access codes with SR can be generalized, both classically and 
quantumly. In particular, one can consider 

• n i-^- 1 codes in base d, d > 2 (called qudits in the quantum case), 

• n m codes with to > 1, 

• n m codes where any k > \ bits (qubits) must be recovered. 

Of course, one can consider several of these generalizations simultaneously. In the setting without 
shared randomness such generalizations have already appeared in the literature (see Sect. 1.2). We will 
briefly introduce the notion of the generalized Bloch vector which we believe can be useful to study 
such generalizations (it has been explicitly used in [5] to prove the impossibility of 2™ m QRAC 
with p> 1/2, when SR is not allowed). 

The notion of the Bloch vector introduced in Sect. 3.1.1 can be generalized for d > 2. For example, 
to write down the density matrix for c? = 3 one uses eight Gell-Mann matrices denoted by Ai instead 
of three Pauli matrices ai defined in equation (35). In general one needs d^ — 1 matrices that span 
the set of all traceless c? x d Hermitian matrices. A convenient choice of A^ are the so called generalized 
Gell-Mann matrices, also known as the generators of the Lie algebra of SU{d), given in [31]. We can 
use them to generalize equation (36): 
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Figure 32: Close-up of the narrow region in Fig. 31 between the quantum upper and lower bound 
(everything is shown relative to the upper bound that corresponds to the horizontal axis). Circles 
indicate the gap between the upper bound and numerical QRACs with SR. Black squares show the gap 
between the upper bound and the lower bound by measurements along coordinate axes (see Fig. 10). 
Dashed line corresponds to the gap between the quantum upper bound and the lower bound by random 
measurements. 
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where A — (Ai, . . . , Xd^-i) and r G M.'^ ^ is the generalized Block vector^^ or coherence vector [28, 30]. 
Since the Ai are chosen so that Tr A^ = and Tr(A,;Aj) = 2Sij, equation (37) generahzes to 

|(Vi|V'2>|' = Tr(pip2) = ^(1 + (d - 1) n • r2). (142) 

To recover a base d digit, we perform a measurement in an orthonormal basis {\ipi) , . . . , |V'd)} of C''. 
Since j'^ = for any pair i ^ j, the corresponding Bloch vectors must satisfy Vi-rj = —jz^- This 

means that they are the vertices of a regular simplex that belongs to a (c? — l)-dimensional subspace 
and is centered at the origin (for d = 2 this is just a line segment). 

On the other hand, in Sect. 4.3 we observed that it might be advantageous to perform measurements 
along orthogonal directions in the Bloch sphere to recover different bits. Let _L Sj be two orthogonal 
Bloch vectors. Then the corresponding quantum states \ilJi) and \(pj) must satisfy \{tfji\ipj)\'^ = ^. This 
is exactly the case when jV'i) and \ipj) belong to different mutually unbiased bases (see Sect. 4.3). 
This suggests that distinct bits should be recovered using mutually unbiased measurements. Note that 
the Bloch vectors of the states from two MUBs correspond to the vertices of two regular simplices in 
mutually orthogonal subspaces. In general, the Bloch vectors of the states from all d + 1 MUBs are the 
vectices of the so called complementarity polytope [32], which is just the octahedron when d — 2. 

The conclusion of Sect. 4.3 and our discussion above suggests the use of MUBs to construct QRACs 
also for d > 2. Such attempts have already been made [6, 8]. Galvao [6] gives an example of 2 i— ^ 1 
QRAC for qutrits {d = 3) and Casaccino et al. [8] numerically investigate (d + 1) i-^ 1 QRACs based 
on MUBs for d-level quantum systems. However, there is a significant difference between the qubit and 
qudit case. Recall that for d — 2 the optimal way to encode the message x is to use a unit vector in 
the direction of (see equation (57) in Sect. 3.4). A similar expression for can be obtained when 
d > 2, but then the matrix p assigned to r = v^/ \\vx\\ according to equation (141) is not necessarily 
positive semidefinite and hence may not be a valid density matrix. However, it is known that for small 
enough values of ||r|| (in our case^^ ||r|| < jrj), all Bloch vectors correspond to valid density matrices 
[29]. Hence, if we cannot use the pure state corresponding to i'^/ ll'i'xll, we can always use the mixed 
state corresponding to -^^zzi'^x/ 11''^ a; II- If one knows more about the shape of the region corresponding 
to valid quantum states, one can make a better choice and use a longer vector, possibly in a slightly 
different direction. Unfortunately, apart from convexity, not much is known about this shape. Already 
for d = 3 it is rather involved [28, 29]. In general the conditions (in terms of the coordinates of the 
generalized Bloch vector r) for p to have non-negative eigenvalues are given in [30, 28]. 

However, for proving only an upper bound, one can ignore all such details. Thus we believe it 
might be possible to generalize our upper bound (see Sect. 3.6 and 3.7) using generalized Bloch vectors. 
It would be interesting to compare such a result with the upper bound (76) that was obtained by 
Ben-Aroya et al. in [16]. 

Finally, another way of generalizing QRACs with SR is to add other resources. A good candidate 
is shared entanglement. 
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A Combinatorial interpretation of sums 

In this appendix we give a combinatorial interpretation of the sums in equations (22) and (23) from 
Sect. 2.3.2. This interpretation is formalized in the form of equations (143) and (144). We referred to 

^^Our normalization follows [30], where the generalized Bloch sphere has radius 1. Another widely used convention is 
to assume radius \j2{d — l)/d, e.g., see [28, 29]. 
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these equations in Sect. 2.3.2 to obtain an exact formula (26) for the average success probabiHty of an 
optimal classical RAC. 

Let us consider a set of n distinct elements and count the number of ways to choose more than half 
of n elements and mark one of them as special. There are two approaches: first choose the elements 
and then mark the special one or first choose the special one and then choose the others. 

In the first scenario there are ways to choose exactly i elements and mark one of them as special. 
If we have to choose more than half, we obtain the sum X]"=m+i ^{T) where m — [^J . 

In the second scenario there are n ways to choose the special element. Then there are I = n — 1 
elements left and at least m of them must be taken to have more than half of n elements in total. The 
number of ways to do this corresponds to the number of subsets of size at least m of a set of I distinct 
elements. Let us consider the cases when / is odd and even separately. 

If n = 2m then I = 2m — 1 is odd. To each "large" subset of size i (m < i < I) we can assign a 
unique "small" subset (the complement set) of size l — i (0 < l — i < m — 1), and vice versa. Each subset 
is either "large" or "small" , so the number of "large" and "small" subsets is the same — it is half of the 
number of all subsets, i.e., 2^/2 = 22™"^ 

If n = 2m + 1 then I — 2m is even. The "large" subsets have m + 1 < * < ^ elements, but the "small" 
ones: 0<l — i<m — 1. Let us call the remaining (^^) subsets of size m "balanced" . In this case the 
bijection between the "large" and "small" subsets holds as well, but it maps the "balanced" subsets 
to themselves. Thus the total number of all subsets is "large" + "small" + (^^) = 2'. The number of 
"large" subsets is (2' + (2^')) /2 = 2^^-' + \ (2^') . 

Both counting methods must give the same results, so for odd and even n we obtain, respectively: 



2m 



2m 



2m- 22^-2 ^ (144) 
We would like to acknowledge Juris Smotrovs for providing this interpretation. 



B POVMs versus orthogonal measurements 

An orthogonal (or von Neumann) measurement is not the most general type of measurement of a 
quantum system. In general a POVM measurement [33, 34] may extract more information. In this 
appendix we show that in the qubit case POVMs can be simulated using a probabilistic combination of 
enhanced orthogonal measurements, as defined in Sect. 3.7 (such a measurement is either an orthogonal 
measurement or a constant function) . To define a POVM we have to introduce the notion of a positive 
semidefinite matrix [35]. 

Definition. A complex square matrix E is called positive semidefinite if ("01 E\il]) >Q for all . 

An equivalent definition is that E is diagonalizable and all eigenvalues of E are real and non-negative. 
Thus E is Hermitian. 

Definition. A positive operator-valued measure (POVM) is a set . . . , E„i\ of positive semidefinite 
matrices such that X]"=i ~ ^- I^^, 34] 

POVM measurements can have an arbitrary number of outcomes, but in the case of n 1 QRACs 
we have to consider only two-outcome single-qubit POVMs. Such a POVM can be specified by {Eq,Ei], 
where Eq is positive semidefinite and Ei — I — Eq. Since Eq is also Hermitian, we can find a basis 
= {\ipo) , IV'i)} in which Eq is diagonal, i.e., Eq — {q°)- In this basis i?i = ( i°f,)- Since both 
Eq and Ei are positive semidefinite, < a < I and < b < I. An arbitrary pure qubit state [i/j) in the 



48 




basis B can be specified by (31). When is measured, the probabihties of the outcomes are 

9 

Po = (^1 Eo \^) = acos^ - + fesin^ -, 

e 9 

Pi = (^1 El = (1 - a) cos^ - + (1 - 6) sin^ -. 

Let us consider the following process (see Fig. 33) that simulates the POVM measurement {Eq, Ei}: 

1. perform an orthogonal measurement in the basis B — {lipo) , 

2. perform the following post-processing of the outcome of the measurement: 

• if the outcome was 0: output with probability a and output 1 with probability I — a, 

• if the outcome was 1 : output with probability b and output 1 with probability 1 — 6. 

To see why this process is equivalent to the POVM measurement {Eq, Ei}, consider a pure qubit state 
given by (31) in the basis B. When j?/;) is measured in the basis B — {|V'o) i IV"!)}; the probabilities 
of the outcomes and 1 are as follows (see also equation (38) in Sect. 3.1.1): 

c 9 

Po = \{lpo\^)f = COS^ 

9 

Pi = KV-ilV")!^ = sin^ 2' 

Now it is simple to verify that the process shown in Fig. 33 has the same outcome probabilities (145) 
as the POVM measurement. However, this process is not a probabilistic combination of enhanced 
orthogonal measurements, since it involves a probabilistic post-processing of the measurement result. 
To obtain the desired result, we have to modify it. The key idea is that with a certain probability the 
output can be produced without performing an actual measurement. 

Let /i = min{a, 6}. Whatever state is input to the process shown in Fig. 33, the probability Pq to 
output is at least /z, because 

Po = Pofl > (po +Pi)m = (147) 

Note that /x does not depend on the state being measured. This means that one can output with 
probability ji without performing an actual measurement. A similar lower bound holds for Pi as well: 

Pi = po(l -a)+ pi{l -b)> (po +Pi)(l - M) = 1 - M, (148) 

where M = max {a, 6} = a + b — fi. Let us consider the following probabilistic combination of four 
decoding strategies: 

• with probability cq: output without performing a measurement. 
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• with prohability ci: output 1 without performing a measurement, 

• with probability cqi: measure in the basis {|V'o) j IV"!)}! 

• with probability ciq: measure in the opposite basis {Itpi) , I'/'o)}- 
The resulting probabihties of outcomes for this process are 

I Po = Co + CoiPo + CioPl , ^^^^^ 
I Pi = Ci + CoiPl + CioPo- 

We can use the lower bounds (147) and (148) for Pq and Pi, respectively, to assign the probabilities cq, 
ci, coi, and cio in the following way: 

= l-(« + ^) + ., (^50) 

= b — fi 

(note that at least one of the probabilities cqi or cio will be zero). It is not hard to verify that after the 
assignment (150) the probabilities Pq and Pi in (149) will match the probabilities of outcomes (145) of 
the POVM measurement. 

Thus for each qubit POVM given by a and b one can find a probabilistic combination of enhanced 
orthogonal measurements given by Cq, Ci, Cqi, and Cio, such that in both cases the probabilities of 
outcomes are the same. 

Example. For a = 6 = 1/2 we have cq = Ci = 1/2 and Cqi = Cio = 0, corresponding to random 
guessing (observe that Eq = Ei in this case). 

Example. However, a = 1 and = corresponds to a projective measurement in basis {\'ipo) , 
because coi — 1 and cio = cq = ci = 0. 

Example. Finally, a — 1 and b — 1 corresponds to a constant function 0, because cq = 1 and 
coi = cio = ci = 0. 
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